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ABSTRACT 


The  purpose  of  this  research  is  to  design  a  digital  phased  array  for  receiving  and 
tracking  a  2.4  GHz  frequency  modulation  (FM)  video  signal  from  an  unmanned  air 
vehicle.  The  tracking  is  done  using  a  monopulse  technique.  Various  numbers  of  elements 
were  simulated  to  access  the  pattern  coverage.  The  beamforming  and  demodulation  are 
performed  digitally  by  a  computer.  Several  antenna  architectures  were  studied  and  the 
requirements  for  hardware  components  such  as  the  analog-to-digital  converters  and 
amplifiers  are  specified.  Several  methods  for  generating  the  carrier  signal  for  the  local 
oscillator  are  also  discussed.  They  include  an  independent  local  oscillator,  a  voltage 
controlled  oscillator  with  feedback  (i.e.,  a  phased  locked  loop)  and  mixing  with  a  phase 
shifted  copy  of  the  received  signal  (Bilotti’s  method). 
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EXECUTIVE  SUMMARY 


As  a  part  of  an  ongoing  project,  this  research  aimed  at  designing  a  digital  phased 
array  for  receiving  and  tracking  a  2.4  GHz  FM  video  signal  from  an  unmanned  air 
vehicle  (UAV).  Tracking  the  signal  from  the  UAV  will  establish  a  data  link  for  data 
transfer  such  as  video,  image  or  audio  between  the  UAV  and  the  ground  station.  The 
tracking  is  done  using  a  monopulse  technique. 

An  array  tradeoff  study  was  done  to  determine  the  benefits  of  1  row  and  2  row 
arrays,  with  8  elements  per  row,  in  terms  of  the  half-power  beamwidths  (HPBWs)  of  the 
azimuth  and  elevation  patterns  of  the  array.  The  one  row  array  was  chosen  primarily  due 
to  its  wide  elevation  pattern  coverage.  This  eliminates  the  need  for  elevation  tracking  as 
the  UAV  is  within  the  HPBW  for  all  UAV  ranges  of  interest.  Adding  a  row  of  8  more 
elements  to  increase  the  signal-to-noise  ratio  (SNR)  will  increase  the  cost  of  equipment. 
The  additional  SNR  is  not  required,  as  the  UAV  transmitter  power  is  sufficiently  high. 
The  sum  and  difference  beams  for  the  Matlab  simulation  of  the  one  row  array  are  verified 
with  the  CST  Microwave  Studio  simulations.  The  printed  circuit  dipole  element  was 
found  to  have  an  E-plane  pattern  HPBW  between  that  of  an  ideal  dipole  and  a  half-wave 
dipole. 

Several  approaches  for  the  demodulation  of  FM  signals  are  explored.  The 
considerations  for  the  choice  of  the  approaches  for  the  FM  demodulation  are  the 
minimum  circuitry  requirement,  the  ability  to  conduct  beamforming,  the  recovery  of  the 
transmitted  signal  and  the  susceptibility  to  phase  errors  and  local  oscillator  (LO)  drift. 
Three  FM  demodulation  techniques  are  considered:  Bilotti’s  quadrature  demodulator,  a 
local  oscillator  with  no  feedback  and  a  digital  phase  locked  loop  (PLL)  quadrature 
detector.  All  three  use  digital  beamforming  and  signal  processing  to  demodulate  the  FM. 
Simulink  models  of  the  three  FM  demodulation  techniques  show  satisfactory 
performance  for  the  recovery  of  a  two-tone  signal. 

The  simulations  were  done  at  a  carrier  frequency  of  50  kHz  instead  of  2.4  GHz,  to 
greatly  save  on  the  processing  time.  A  LO  with  no  feedback  was  chosen,  as  it  uses  the 

xiii 


least  circuitry  with  satisfactory  performance  for  the  recovery  of  the  two-tone  signal. 
Beamforming  can  also  be  performed  easily  as  the  raw  in  phase  (I)  and  quadrature  (Q) 
data  are  preserved.  This  technique  also  eases  the  computer  processing.  Since  no  error 
correction  needs  to  be  computed  for  the  LO  as  it  does  for  a  PLL.  Simulations  were  also 
done  with  relative  phase  errors  introduced  in  the  channels  and  a  frequency  offset  in  the 
LO  without  feedback.  The  Simulink  simulations  indicate  that  the  recovery  of  the  two- 
tone  signal  can  still  be  achieved  with  the  expected  errors. 

Finally,  the  possible  hardware  components  are  briefly  introduced  and  the 
minimum  power  for  the  receiver  module  was  used  to  determine  the  maximum  operating 
range  of  the  array  using  the  one-way  Friis  equation. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Phased  array  development  began  around  1958,  when  early  radio  transmitters  and 
radars  used  multiple  radiating  elements  to  achieve  the  desired  antenna  radiation  patterns 
[1].  As  technologies  advance,  existing  systems  are  improved  and  upgraded.  Sometimes, 
new  systems  are  created  when  the  existing  systems  become  obsolete.  Unmanned  systems 
and  array  antennas,  two  technologies  that  improve  continuously,  have  brought  many  new 
capabilities  to  military  operations.  The  superiorities  of  array  antennas  are  the  digital  and 
computerized  processing  capabilities  that  can  be  built  into  them.  Digital  phase  shifting  to 
form  a  beam,  the  electronic  and  rapid  beam  scanning  ability,  and  its  broadband  potential 
are  the  advantages  of  array  antennas. 

Since  the  Gulf  War  in  the  1990s,  the  military  has  a  new  awareness  of  the 
importance  of  technological  advances  such  as  stealth  ability,  unmanned  systems  and 
global  positioning  systems  (GPS)  in  improving  the  success  of  a  war.  One  of  these 
technologies  is  the  digital  phased  array.  It  has  many  advantages,  such  as  scanning  rapidly 
without  physically  moving.  Multiple,  simultaneous  beams  can  be  achieved  by  the  means 
of  digital-beam  forming  (DBF).  Error  correction  and  self-calibration  can  be  achieved 
easily  in  the  signal  processing.  DBF  also  allows  pattern  control  to  direct  the  nulls  in  the 
pattern  in  the  direction  of  noise  jammers,  and  also  suppresses  clutter,  which  is  a  cause  of 
degradation  in  range.  DBF  also  gives  the  opportunity  to  compensate  for  the  failure  of  the 
elements  in  the  array.  Different  data  rates,  depending  on  the  operating  environment,  can 
be  achieved  with  DBF.  The  array  antennas  can  be  used  in  many  applications,  such  as 
communications,  data-links,  radar  (search  and  track),  and  electronic  warfare  (EW).  This 
shows  perhaps  the  most  important  advantage  of  these  antennas:  the  multi-functionality. 

Unmanned  aerial  vehicles  (UAVs),  one  important  type  of  unmanned  system,  play 

a  very  critical  role  in  military  operations.  UAVs  have  a  variety  of  uses,  such  as 

surveillance  and  reconnaissance,  target  recognition,  battle  damage  assessment,  EW, 

search  and  rescue,  and  traffic  monitoring.  All  of  these  can  be  conducted  without 

1 


endangering  the  life  of  a  pilot,  as  they  do  for  manned  systems.  Important  roles  for  the 
UAVs  are  troop  support  and  carrying  out  reconnaissance  and  surveillance  missions, 
which  require  maintaining  a  data-link  with  troops  in  order  to  send  any  data  collected, 
such  as  video  images  or  audio,  back  to  the  command  center. 

During  operations,  it  is  imperative  to  continuously  maintain  a  data  and  control 
link  with  the  operator.  This  requires  the  ground  station  antenna  to  track  the  UAV  so  the 
antenna  beam  is  pointed  properly.  The  purpose  of  this  research  is  to  design  and  build  an 
array  to  angle-track  a  UAV  and,  eventually,  to  accomplish  the  data  transfer  from  the 
UAV  to  the  ground  station.  The  antenna  type  on  the  UAV  is  not  a  subject  of  this  research 
and  it  will  be  assumed  to  be  a  low-gain  omni-directional  antenna.  What  is  of  concern  is 
the  array  system  that  will  be  designed  to  achieve  the  purpose  that  was  just  mentioned. 

Why  is  an  array  antenna  selected  for  this  application?  The  narrow  beam  of  the 
antenna  must  be  kept  on  the  UAV  as  it  moves  and  maneuvers.  The  beams  of  an  array  can 
track  without  physically  moving  the  antenna.  Digital  array  systems  can  be  reconfigured 
and  different  beamwidths  can  be  used,  depending  on  the  UAV  path  and  the  data  rate 
desired  for  the  type  of  operation.  Array  antenna  systems  can  track  and  follow  more  than 
one  UAV  at  a  time  (multiple  tracking)  [2],  Figure  1  shows  a  High  Gain  Advanced  GPS 
Receiver  (HAGR)  phased  array  antenna. 


Figure  1.  HAGR  16  element  antenna  array.  From  [3]. 
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B. 


PHASED  ARRAY  TECHNOLOGY  FOR  MILITARY  APPLICATIONS 


Since  1901,  array  antenna  technology  has  made  tremendous  progress  with  types 
that  vary  from  linear  to  planar  to  conformal.  Digital  array  applications  enable  the  main 
beam  of  arrays  to  be  steered  electronically,  rather  than  mechanically.  A  phased  array  has 
its  radiating  elements  excited  by  the  relatively  phased  currents  and  its  main  beam  can 
scan  the  desired  field  continuously.  Since  the  1980s,  phased  arrays  have  gained  wider  use 
in  radar  and  communication  applications. 

Phased  array  radars  are  used  in  a  number  of  navy  ships,  such  as  the  US  Navy 
Aegis  ships,  the  German  Navy’s  FI 24  frigates  and  LCF  frigates  of  the  Royal  Netherlands 
Navy  [4],  The  advantages  of  the  phased  array  radars  used  in  communication  systems  are 
summarized  below: 

1.  Electronically  scanning  the  beam  and  its  multiple  element  system 
increases  the  reliability  of  the  radar  system,  and  maintenance  is  simplified  by 
replacing  the  faulty  elements  with  new  ones  without  a  need  to  return  to  the 
shipyard.  An  antenna  can  have  a  significant  number  of  its  elements  destroyed 
and  still  be  able  to  function,  albeit  at  a  reduced  capability.  Electronic  arrays 
have  no  moving  parts,  do  not  suffer  from  vibration  and  are  inherently  better 
stabilized  against  ship  or  aircraft  motion,  as  the  electronic  array  simply  steers 
the  main  beam  counter  to  the  detected  motion. 

2.  It  is  suitable  for  low  observable  platforms  as  antenna  elements  of  the  array, 
like  a  microstrip  patch  antenna,  can  conform  to  the  structure  of  stealthy 
platforms.  Its  agile  scanning  rate,  with  a  narrow  beam  and  lower  sidelobes, 
gives  a  low  probability  of  interception  and  is  less  vulnerable  to  jamming. 

3.  The  integrated  circuit  (IC)  technology  in  microwave  monolithic  ICs 
(MMICs)  is  being  driven  by  the  need  to  integrate  disparate  technology 
functions  into  one  unit.  Through  the  use  of  IC  technology,  phased  arrays  can 
reach  improved  cost,  size,  and  weight  metrics  that  have  previously  prevented 
the  widespread  adoption  in  mass  market  scenarios.  The  high  power  handling 
capability  and  low  loss  of  ferrite  phase  shifters  allow  power  amplifiers,  low 
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noise  amplifiers  and  phase  shifters  to  be  chosen  to  achieve  high  processing 
performance,  with  relatively  little  impact  to  cost.  Hence  the  solid  state  and  IC 
technologies  will  benefit  the  phased  array  in  terms  of  costs,  weight  and 
reliability  [5,  6], 

4.  Its  ability  to  create  multifunction  systems  (radar,  sensing,  communications, 
etc.)  with  a  single  device  enabled  by  the  system’s  software-driven  nature 
increases  the  tactical  flexibility  of  its  operational  employment.  This  capability 
is  explained  in  the  following  scenarios: 

•  The  radiated  power  from  each  element  can  be  adjusted  when  each 
element  has  its  own  transmitter.  At  the  expense  of  cost,  the  controlling 
of  the  radiated  power  in  sectors  addresses  the  concern  for  emission 
control  and  co-site  interference. 

•  The  elements  can  be  programmed  to  communicate  with  multiple 
UAVs  by  scanning  one  single  beam  rapidly  or  producing  multiple, 
independent  beams.  This  would  increase  the  dwell  time  on  each  UAV. 

•  The  multifunction  system  reduces  the  number  of  antennas  on  the 
military  platforms,  hence  reducing  weight  and  the  radar  cross  section 
of  the  platform. 

•  Its  flexibility  in  electronic  control  allows  for  customized  software 
solutions  to  address  different  user  scenarios  with  one  platform  [5], 

•  A  graceful  degradation  is  allowed  when  some  of  the  elements  are 
damaged  enabling  the  array  system  to  continue  functioning. 

The  frequency  bandwidth  is  a  major  design  challenge  for  the  array  system  as 
bandwidth  is  a  scarce  commodity  in  modem  communications  systems,  with  the  Federal 
Communications  Commission  (FCC)  deciding  which  parts  of  the  spectrum  can  be  used 
for  emerging  communication  standards.  Despite  the  multi-functionality  of  the  current 
arrays,  different  functions  need  different  bands,  due  to  the  unique  characteristics  and 
behaviors  of  the  electromagnetic  waves  in  different  frequency  bands.  Multifunction 
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arrays  can  be  used  for  surveillance,  tracking,  communications,  direction  finding  and 
weapon  guidance  purposes.  The  radar  should  be  designed  in  lower  bands  for  longer  range. 
The  AN/SPY- ID  operates  in  S  band.  Current  electronic  support  equipment  operates 
within  the  range  from  2  to  20  GHz.  However,  due  to  the  dense  electronic  environment,  it 
is  very  difficult  for  the  electronic  warfare  operators  to  detect  radars  operating  in  the  UHF 
band  in  the  littoral  waters. 

The  complexity  of  the  phased  arrays  will  increase  the  cost.  Commercial  off-the- 
shelf  (COTS)  hardware  is  to  be  exploited  to  reduce  the  cost  of  the  phased  arrays.  The 
ultimate  aim  of  the  antenna  engineer  is  to  decrease  the  cost  of  the  phased  array  antennas 
without  sacrificing  performance  [2], 

C.  PREVIOUS  WORK 

This  research  is  part  of  the  ongoing  project  that  started  with  the  thesis  by  Bartee 
[7].  Esswein  [8]  and  Gezer  [2]  continued  the  work  and  demonstrated  the  feasibility  of 
building  a  digital  phased  array  with  available  COTS  hardware. 

Bartee  designed  and  constructed  the  distributed  phased  array  at  X-band. 

Esswein  constructed  a  24-element,  three-dimensional,  dual-ground  plane,  active 
phased  array.  A  Genetic  Algorithm  (GA)  was  used  to  determine  the  element  excitations 
and  the  measured  radiation  beam  was  in  agreement  with  the  theoretical  calculations.  The 
modulator  feeding  each  printed  antenna  element  used  for  this  array  was  the  Analog 
Devices  AD8346EVAL  quadrature  modulator  board,  which  is  now  being  incorporated 
into  a  transmit/receive  (T/R)  module  for  radar  applications. 

Gezer  constructed  a  1  by  8  phased  array  to  track  a  signal  transmitted  continuously 
by  a  UAV  operating  at  2.4  GHz.  This  array  used  a  monopulse  difference  beam  to  track 
the  UAV  signal.  The  array  was  designed  based  on  typical  UAV  operational  requirements, 
which  led  to  the  selection  of  the  array  parameters,  such  as  the  number  of  elements  and 
spacing. 
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D.  SCOPE  OF  RESEARCH 


State  of  the  art  technology  in  electronic  systems  is  a  means  to  increase  the 
information  superiority  and  provide  the  decision  makers  with  a  well  described  tactical 
picture  of  the  war  scene.  UAVs  can  carry  out  many  reconnaissance  and  surveillance 
missions.  They  also  provide  the  necessary  information  to  the  troops  by  being  integrated 
into  the  network-centric  warfare  elements  as  seen  in  Figure  2. 


Data  Link  Relay 


Figure  2.  The  cooperation  among  the  elements  of  a  UAV  system.  From  [10]. 

As  an  element  of  the  network,  the  ground  control  stations  of  UAV  systems  can 
exploit  the  phased  array  technology.  Array  antennas  are  a  very  sophisticated  means  to 
accomplish  accurate  UAV  tracking  and  maintain  a  good  data-link. 

The  research  goal  was  to  design  a  digital  phased  array  for  tracking  and  receiving 
FM  video  signals  transmitted  at  2.4  GHz  from  a  UAV.  The  design  incorporates  COTS 
hardware  to  reduce  cost.  Simulation  is  another  way  to  reduce  the  cost  of  the  design  and 
comprehend  the  consequences  of  design  changes  before  manufacturing,  and  to  see  if 
there  needs  to  be  any  adjustment.  Thus,  the  design  of  the  tracking  array  antenna  system  is 
verified  through  simulations  before  being  constructed  [2], 
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E.  ORGANIZATION  OF  THESIS 


Chapter  II  discusses  the  concepts  of  arrays.  The  concepts  of  a  digital  antenna  and 
digital  beamforming  are  explained,  and  the  mathematical  theory  for  arrays  with  ground 
planes  is  elaborated.  Amplitude  comparison  monopulse  systems,  and  how  tracking  of  the 
UAV  is  achieved,  along  with  external  causes  of  errors,  are  mentioned. 

Chapter  III  addresses  the  results  and  analysis  of  the  array  tradeoff  study.  It 
discusses  the  results  obtained  for  the  1  row  and  2  row  arrays  with  8  elements.  Also,  the 
effect  of  polarization  is  studied. 

Chapter  IV  addresses  the  FM  demodulation  techniques  being  considered.  Several 
methods  for  generating  the  carrier  signal  are  also  discussed.  They  include  an  independent 
local  oscillator,  a  voltage  controlled  oscillator  with  feedback  (i.e.,  a  phased  locked  loop) 
and  mixing  with  a  phase  shifted  copy  of  the  received  signal  (Bilotti’s  method). 

Chapter  V  discusses  the  system  architecture  of  the  array  and  its  considerations. 
Hardware  implementation  using  COTS  components  is  also  described. 

Chapter  VI  summarizes  what  was  done  in  this  research  and  gives 
recommendations  for  future  work. 
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II.  TRACKING  ARRAY  DESIGN 


A.  WHY  AN  ARRAY  ANTENNA? 

An  antenna  array  is  a  configuration  of  individual  radiating  elements  that  are 
arranged  in  space  to  produce  a  directional  radiation  pattern.  Single-element  antennas  have 
broad  and  low  directivity  radiation  patterns,  which  are  not  suitable  for  long  distance 
communication  or  accurate  tracking.  High  directivity  can  be  achieved  by  increasing  the 
electrical  dimensions  of  the  single-element  antennas  in  terms  of  wavelength,  hence  its 
physical  size.  This  makes  the  antenna  too  cumbersome  for  mobile  operations.  The 
alternative  is  to  combine  multiple  electrically  small  antennas  to  form  an  array.  Antenna 
arrays  come  in  various  geometrical  configurations,  including  linear,  planar  and 
volumetric.  The  planar  antenna  array  is  often  used  in  communication  and  radar 
applications,  and  therefore  this  general  configuration  will  be  analyzed  [11]. 

Conventional  phased  arrays  employ  microwave  circuits  to  distribute  signals  to  the 
elements  or,  on  receive,  collect  signals  from  the  elements.  These  beamforming  networks 
(or  beamformers),  as  shown  in  Figure  3,  can  be  bulky,  narrowband  and  lossy.  Some  of 
these  problems  can  be  avoided  by  using  digital  beamforming. 


Figure  3.  Conventional  microwave  beamforming  network.  From  [12], 
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B.  DIGITAL  ARRAYS1 


The  basic  architecture  of  a  digital  array  is  shown  in  Figure  4.  For  a  two-way 
communication  application,  each  element  must  be  capable  of  transmission  and  reception. 
In  the  transmit  channel,  the  waveform  is  generated  and  amplified  before  passing  to  the 
array  element.  On  receive,  the  signal  is  amplified  and  detected.  The  detected  signals  are 
collected  by  a  processor  which  forms  the  beam  and  performs  other  communication 
functions. 

In  order  for  all  of  the  elements  to  send  and  collect  data  coherently,  a  common 
frequency  (phase)  reference  must  be  distributed  to  all  of  the  elements.  This  signal  is 
referred  to  as  the  local  oscillator  (LO)  signal.  The  elements  must  also  be  time 
synchronized  so  that  transmission  can  occur  simultaneously  and  the  transmitted  signals 
arrive  at  the  other  end  of  the  link  at  the  same  time. 


LO  and  Timing  Tx  channel:  Waveform  Generator, 


Generally  the  array  has  a  modular  architecture  where  the  transmit  (Tx)  and 
receive  (Rx)  functions  are  channeled.  The  transmit  portion  generates  the  waveform, 

1  Much  of  the  material  in  this  section  is  adapted  from  unpublished  notes  by  Prof.  D.  Jenn. 
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amplifies  and  filters  it  before  passing  it  to  the  antenna.  In  modem  arrays,  the  waveform  is 
generated  using  a  direct  digital  synthesizer  (DDS).  In  its  simplest  form,  a  DDS  is  a  clock 
and  memory  with  a  stored  waveform.  The  waveform  is  read  out  at  a  rate  determined  by 
the  clock.  The  DDS  output  is  usually  a  baseband  signal  that  has  to  be  upconverted, 
although  recent  DDS  clock  speeds  are  achieving  2  GHz.  Therefore,  at  VHF  radar 
frequencies  and  below,  it  may  be  possible  to  amplify  the  DDS  output  directly  without 
upconversion.  For  the  tracking  application,  this  research  is  concerned  only  with  the 
receive  function.  Therefore,  the  remaining  discussion  only  addresses  the  receive  channel. 

On  the  receive  side  the  returned  signal  is  amplified  with  a  low  noise  amplifier 
(LNA),  downconverted,  filtered  and  then  sampled.  State  of  the  art  analog-to-digital 
converters  (ADCs)  are  in  the  2  GHz  range,  so  at  VHF  and  below  it  may  be  possible  to 
sample  the  waveform  directly  without  downconverting. 

The  in-phase  (7)  and  quadrature  (Q)  form  of  the  data  are  handled  by  the  hardware. 
If  the  signal  at  carrier  frequency  coc  =  In  fc  is: 


s(t )  =  A(t )  cos  [2  n  fct  +  <E>(t)] 


(1) 


the  in-phase  and  quadrature  representation  is: 
s(t)  =  I{t)  cos(2nfct)  -  Qit)  sini2nfct) 

where 

lit)  =  Ait)  cos  [0(0] 

8(0  =  A(t)  sin  [0(0] 


Ait)  = 
®it)  = 


V i  if )2 +£?o 
Q(0 


tan 


-1 


1(0 


(2) 


(3) 


(4) 


The  signal  is  conveniently  shown  in  the  HQ  plane  (Figure  5),  where  the  distance  from  the 
origin  is  Ait)  and  the  rotation  from  the  /  axis  is  <f>(t) . 
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Figure  5  also  shows  a  quadrature  demodulator,  the  Analog  Devices  AD8347  [14], 
which  operates  from  800  MHz  to  2.6  GHz  with  a  60  MHz  instantaneous  bandwidth.  The 

demodulator  has  differential  outputs  (7+,  I~,  Q+,  and  Q~  )  for  obtaining  the  baseband 
in-phase  and  quadrature  components  of  the  received  signal.  The  components  are  given 
by: 


=  I+  -  I 

=  Q+-Q~ 


The  hardware  aspects  are  discussed  further  in  Chapter  V,  Section  A. 


Figure  5.  HQ  signal  diagram  (left)  and  the  AD8347  quadrature  demodulator  (right). 

From  [12]. 


The  demodulation  process  is  analyzed  using  the  complex  envelop  of  the  signal 
u(t)  [15]: 

u(t)  =  I(t)  +  jQ(t) 

s(t )  =  Re  jw(t)e^c?  j  (6) 

=  ^ \u(t)ej0)ct  +  u  *  (t)e~jc°ct  ]  • 
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The  mixing  process  (shown  in  Figure  6)  is  equivalent  to  multiplying  the  received  signal 
with 

COS  COLO  =  ^(ej<°Lot  +  e-j®LOt  )  (7) 

in  the  /  channel  and  with 


-sm  o)LO  = 
in  the  Q  channel. 


—  _e~JC0LOt ) 


(8) 


KO 


Figure  6.  Quadrature  demodulation  process.  From  [12]. 


For  the  in-phase  channel: 

rit )  cos  t oLO  =  ^  [u(t)eJ0^  +  u  *  (t)e~J6^  ]  {e’0^  +  ) 

=  u(t)  cK(»r+6>i.n)t  !  J(ojr-coin)t 

4  4 

(9) 

I  u*(0  c-/(a>r-am)t  ,  u*(t)  C-J'(c»r+C0j.n)t 

4  4 

=  Re  {*)  e7(®c  +®io )?  |  +  Re  e-/(®c  -®lo )*  | 

where  Re{.}  is  the  real  operator  and  the  *  is  complex  conjugate. 

If  direct  down  conversion  is  used  (no  intermediate  frequency,  coLO  =  coc ),  then  the 
first  term  is  at  frequency  2 a>c ,  which  is  filtered  by  the  low  pass  filter  (LPF).  The  second 
term  is  the  desired  base  band  signal. 
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A  similar  analysis  for  the  quadrature  channel  gives  two  terms: 

-r(0sin®£O  =im|^eA®c-®ioV|_im|^e-/(®c+®LO)^|-  (10) 

where  Im{.}  is  the  imaginary  operator  and,  similar  to  the  I  channel,  the  second  term  is  at 
frequency  2 coc  and  is  filtered.  The  /  and  Q  outputs  are  combined  in  the  processor  to 
form  the  complex  signal  [12]. 

C.  DIGITAL  BEAMFORMING 

For  a  simple  linear  array  along  the  z  axis  in  Figure  7,  the  array  factor  is  given  by: 

N  N 

AF{6,  <f>)  =  X  >'/;  exP  \j fid n  cos  #]  =  X  AneJCCn  exp  [jj3dn  cos  <9]  (11) 

n=  1  n= 1 


Figure  7.  Digital  array  beamforming  on  receive.  From  [12]. 

where  N  is  the  number  of  elements,  J3  =  In / A  (A  is  wavelength),  wn  =  Aneja'n  is  a 
complex  weight  added  by  the  processor  and 

[2»-(JV  +  l)] 

un  2  u 

is  the  location  of  element  n  relative  to  the  origin  at  the  center  of  the  array. 
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(12) 


The  demodulators  provide  spatial  samples  of  the  exponential  factor  in  Equation 


(11): 


In  =  cos  [  j /3dn  cos  <9] 

(13) 

Qn  =  sin  [./Mr  COS  #] 

(14) 

Therefore,  to  obtain  a  response  equivalent  to 

a  beam  scanned  in  the  direction  6S 

weights  must  be  given  by: 

wn  =  exP  [~jPdn  COS  0S  ] 

(15) 

Note  that  compensation  for  any  phase  errors  in  the  hardware  could  be  included  in  the 
exponent  of  wn  .  They  can  be  obtained  by  calibration  [12]. 

D.  ARRAYS  WITH  GROUND  PLANES 

This  section  derives  the  formulas  used  in  the  Matlab  codes  to  obtain  the  radiation 
pattern  of  a  two-dimensional  array  with  a  rectangular  element  lattice.2  Figure  8  shows  a 
planar  array  of  point  sources  in  the  x-z  plane. 


Figure  8.  Planar  array  centered  at  the  origin.  After  [16]. 


2  This  analysis  is  updated  from  unpublished  notes  by  Prof.  David  C.  Jenn  of  NPS. 
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The  pattern  angles  0  and  <f)  and  the  array  orientation  are  shown  in  Figure  9. 


Figure  9.  Array  orientation  and  pattern  angles  6  and  </> .  After  [17]. 


If  the  array  is  centered  at  the  origin,  then  the  element  locations  can  be  written  as: 
2n  -  {Nx  + 1) 


xn  = 


zm 


2 

2m  -  (Nz  + 1) 


dx,  fi  l,...,Nx 


(16) 


dz,  m  =  \,...,Nz 


where  dx  and  dz  are  the  lattice  dimensions  and  Nx  and  Nz  are  the  number  of  elements 
along  x  and  z,  respectively.  The  array  factor  is: 


AF  = 


A  JPxn  sin  6  cos  <l!jPzm  cos  9 
^ mn e  e 


n- 1  m- 1 


(17) 


If  the  feeding  arrangement  results  in  a  separable  distribution  (which  is  common) 
then  Amn  =  AmAn  and 


AF  =  AFXAFZ 


Nr 


^  AneJPx"  sin9cos $  Ame^z,n  cos  9 


n-\ 


m- 1 


(18) 
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Furthermore,  if  phases  are  introduced  to  scan  the  beam  to  the  direction  ( 6S,(/>S ) 


then 


AF  = 


Nz 

AneJPXn  (sin  0 cos ^_sin  9s  cos  & } 


A  JPzm  (COS  0-COS  <9V) 


n= 1 


(19) 


Finally,  if  the  array  is  uniformly  excited  ( Amn  =  1 ),  the  sum  becomes  a  geometric  series 
with  the  following  closed  form  result  [18]: 


AF  =  AFXAFZ 


sin(A^x  /2)  sin(jVz^z  /2) 

sin(^/2)  sin(^z/2) 


(20) 


where 


y/x  =  pdx  (sin  6  cos  (j)  -  sin  6S  cos  <j>s ) 
i//z  =  fid.  (cos  6  -  cos  0S ) 

Now,  if  the  array  is  placed  at  height  h  above  a  perfect  electric  ground  plane  of 
infinite  extent,  then  the  method  of  images  can  be  used  to  compute  the  field  above  the 
ground  plane.  The  total  array  (sources  plus  images)  can  be  considered  as  a  two-element 
array  along  they  axis,  where  each  element  itself  is  a  planar  array  as  shown  in  Figure  10. 


Figure  10.  Two  element  linear  array  along  they  axis,  where  each  element  is  a  planar 

array.  After  [16]. 
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The  ground  plane  factor  (GF)  for  a  two-element  linear  array  along  the  y  axis,  with 
the  images  out  of  phase  with  the  sources  is: 

GF  =  ejPhsmdsin<fi  _e-j/3hsin0sm<fi  =  2 jsin(0hsin0sinfl  (21) 

The  total  normalized  pattern  factor  is  obtained  using  the  principle  of  pattern 
multiplication 


F(8,<P)  =  (AFX-AFZ) 

norni '  ^^norm  '  ^^nonn 


=  (AFxAF7GF-EF ) 

v  x  z  /norm 

_  sin(A^x  /  2)  sin(  Az^z  /  2) 

Nx  sin  (Vx  /  2)  Nz  sin  (Vz  1  2) 


sin  (j3h  sin  6  sin  <f>)  EFnorm 


where  EFnorm  is  the  normalized  element  factor. 


(22) 


In  general,  the  element  factor  will  have  both  6  and  (f>  components,  and  therefore 
F  has  both  components  as  well.  Referring  to  Figure  9,  the  dipoles  are  aligned  with  the  x 
axis  (vertical  polarization).  For  a  half-wave  dipole  along  the  x  axis,  the  element  factors 
are  (see  Appendix  A): 


EF/i  = 


Flt  Jm  e 


~iPr 


2  nr 


EF^=- 


~iPr 


2  nr 


cos|ySin^cos^ 
1-sin2  #cos2  (j) 
coslysin^cos^ 
1-sin2  ^cos2  (j> 


cos  6  cos  (j) 


sin  <j> 


(23) 


To  normalize  these  equations  the  leading  factor 


Mm  e~Jflr 

2  7tr 


is  simply  removed. 


If  the  dipoles  are  ideal  (with  constant  current  and  length,  /  «  1),  instead  of  half 
wave  dipoles  then  the  terms  in  square  brackets  [  ]  are  reduced  to  1  and  /  is  added  to  the 
leading  factor.  Therefore,  after  normalization  of  the  element  factors,  only  the 
trigonometric  functions  to  the  right  of  the  bracketed  terms  remain  [14].  This  concept  is 
used  in  the  Matlab  codes  as  shown  in  Appendix  B. 
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E.  AMPLITUDE  COMPARISON  MONOPULSE  SYSTEM 


The  tracking  technique  for  the  UAV  application  is  an  amplitude  comparison 
monopulse  system  [2],  Amplitude  monopulse  direction  finding  is  a  common  technique 
employed  in  radar  systems  to  improve  the  accuracy  with  which  the  direction  of  the  arrival 
of  a  single  pulse  can  be  estimated.  Two  “squinted”  beams  are  formed  simultaneously, 
such  that  they  overlap  at  the  3  dB  points  of  the  beams,  as  shown  in  Figure  11. 


Figure  1 1 .  Two  overlapping  beams.  From  [20], 

By  comparing  the  relative  amplitude  of  the  voltages  of  the  signal  being  tracked 
out  of  the  two  beams,  its  position  in  the  beams  can  be  determined  with  an  accuracy 
dependent  on  the  signal- to-noise  ratio  (SNR).  A  common  rule  of  thumb  is  that  the 
accuracy  of  a  tenth  of  beamwidth  can  be  achieved  with  an  SNR  of  10  dB  [18],  The 
difference  of  the  amplitudes  obtained  from  these  two  beams  (the  difference  beam)  gives 
the  angular  error.  If  the  phase  of  the  sum  pattern  is  compared  to  that  of  the  difference 
pattern,  the  direction  of  the  angular  error  is  found,  as  illustrated  in  Figure  12. 


z 

Antenna  (top  view) 

-A(Az)  ^ 

erro^ 

+ A(Az)^r 

Figure  12.  Illustration  of  azimuth  sum  and  difference  beams.  After  [20], 
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The  sum  beam  is  used  for  transmission  and  the  difference  and  sum  beams  are 
used  together  for  the  receive  side  in  monopulse  systems.  The  sum  and  difference  beams 
are  also  shown  in  Figure  13. 


-15  -10  -5  0  5  10  15 

Bearing 


Figure  13.  Sum  and  difference  beam  output  voltages.  From  [20], 


The  ratio  of  the  difference  beam  to  the  sum  beam  is  used  to  generate  the  angle 
error  difference.  The  tracking  process  uses  the  complex  voltage  ratio: 

A  _  difference  voltage  /0/n 

2  sum  voltage 


The  complex  voltage  in  Equation  (25)  can  be  separated  into  the  real  and 
imaginary  terms: 

Rej|j  =  |dCOS<5'  (25) 


Im{f}  =  j^fSin<^  (26) 

where  8  is  the  relative  phase  between  sum  and  difference  channels.  Only  the  Re  j-^-j  is 

processed,  as  it  defines  the  side  of  the  null  with  the  sign  of  the  ratio  and  the  signal 
information  exists  only  in  the  real  term  of  the  ratio  [17]. 


20 


The  plot  of  the  pattern  angle  versus  normalized  difference  to  sum  voltage  ratio  is 
shown  in  Figure  14. 


Figure  14.  Normalized  difference  to  sum  voltage  ratio  versus  the  pattern  angle. 

From  [20]. 

The  slope  in  the  linear  region  creates  an  error  signal  that  helps  the  difference  null 
to  correct  the  beam,  such  that  the  target  signal  is  within  the  null.  This  ratio  is  proportional 
to  the  multiplication  of  the  slope  constant  ks  with  the  pattern  angle  in  degrees  as: 

f  (27) 

where  0B  is  the  half-power  beamwidth  (HPBW).  The  principle  mentioned  above  is 

exploited  in  the  tracking  system  of  this  array  application.  It  must  also  be  noted  that  the 
accuracy  of  this  techniques  can  be  affected  by  unwanted  multipath  reflections  [2], 

F  OTHER  EXTERNAL  CAUSES  OF  ERROR 
1.  Multipath 

This  phenomenon  occurs  when  a  signal  is  reflected  from  the  external 
environment,  such  as  the  terrain  and  other  objects.  The  signal  reaches  the  receiver  from 
different  paths  and  enters  the  antenna  from  different  points  on  its  main  beam  or  sidelobes. 

The  elevation  tracking  errors  are  contributed  to  by  smooth  highly  reflecting  ground 
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surfaces.  Diffuse  scattering  from  buildings  or  rough  surfaces  cause  both  azimuth  and 
elevation  tracking  errors.  The  multipath  errors  increase  as  the  altitude  of  the  target  gets 
lower. 


A  rms  value  of  the  multipath  error  is  calculated  from  [21]: 


°E  = 


P@B 

yJ&Gge  ( peak ) 


(28) 


where  p  is  the  magnitude  of  the  surface  reflection  coefficient,  0B  is  antenna  half-power 
beamwidth  (HPBW)  and  Gse(peak )  is  the  power  ratio  of  the  tracking  antenna  sum- 

pattern  peak  to  the  error-pattern  peak  sidelobe  level  at  the  angle  of  arrival  of  the  image 
signal.  One  of  the  objectives  in  looking  at  a  planar  array  is  to  reduce  the  effect  of 
multipath.  Adding  more  rows  decreases  the  elevation  HPBW  ( 0B )  and  also  decreases 

c 7E  according  to  Equation  (29). 


2.  Crosstalk  Caused  by  Cross-Polarized  Energy 

Crosstalk,  or  cross-polarization  coupling,  is  due  to  the  cross-polarized  energy  of 
the  source,  like  the  undesired  capacitive,  inductive  or  conductive  coupling  from  one 
circuit,  part  of  a  circuit  or  channel  to  another.  This  effect  is  generally  negligible  with  a 
good  array  antenna  element  design.  However,  in  situations  where  crosstalk  is  significant, 
a  large  tracking  error  will  result  in  the  loss  of  track. 

3.  Troposphere  Propagation 

In  the  troposphere,  which  includes  altitudes  between  0  through  10  km  of  the 
surface,  the  variation  of  the  air  density,  dielectric  constants  and  the  thickness  of  the 
clouds  causes  beam  bending.  Tropospheric  refraction  can  lead  to  a  difference  in  the 
sensed  direction  from  the  bent  ray  and  the  actual  direction.  At  standard  atmosphere 
conditions  and  short  range,  this  error  is  not  significant  [2], 
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G.  SUMMARY 


Normally  in  a  digital  array,  each  element  can  act  as  a  transmitter  and  receiver.  A 
LO  is  used  to  synchronize  the  elements  to  enable  concurrent  transmission  and  reception 
of  signals.  The  digital  beamforming  concept  is  discussed.  The  ground  plane  effect  on  the 
radiation  pattern  due  to  the  image  created  by  the  ground  plane  is  derived.  The  tracking 
technique  used  is  an  amplitude  comparison  monopulse  system  where  the  sum  and 
difference  beams  help  in  accurately  tracking  the  UAV.  In  the  next  chapter,  the  computed 
array  patterns  are  compared  to  simulated  pattern  using  CST  Microwave  Studio. 
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III.  ARRAY  TRADEOFF  STUDY 


This  chapter  addresses  the  antenna  requirements  for  the  array  design.  The  antenna 
requirements  are  the  scan  range,  HPBW  and  coverage  in  elevation  based  on  the  UAV 
range  and  altitude. 

A.  HALF-POWER  BEAMWIDTHS 

Several  array  antenna  analyses  were  done  to  determine  the  suitable  array 
configuration  for  the  tracking  function.  Gezer  [2]  designed  and  built  a  linear  array  of  8 
elements.  The  question  was  whether  increasing  the  number  of  elements  to  form  a  planar 
array  would  improve  the  tracking  and  range  performance.  First,  a  2  by  8  element  array  is 
compared  to  the  original  1  by  8  element  array  to  determine  the  worthiness  of  using  a 
larger  number  of  elements.  Next,  the  difference  and  sum  beam  patterns  are  compared 
with  the  results  generated  using  the  CST  Microwave  Studio  simulations  from  the  thesis 
by  Gezer. 

Simulations  used  the  following  parameters:  Nx  =  1  or  2,  Nz  =  8,  dx=  10  cm, 
dz  =6.5  cm,  h  =  0.25 X  .  The  azimuth  and  elevation  patterns  were  computed  for  2.4  GHz 
using  the  Matlab  codes  in  Appendix  B.  Figure  15  shows  the  elevation  patterns  of  arrays 
with  1  and  2  rows  ( Nx  =  1 ,  2).  The  half-power  angles  are  indicated  by  the  arrows  with 

the  degrees  shown.  From  Figure  15,  the  half-power  beamwidth  of  a  1  row  array  is 
80  °  and  the  HPBW  of  a  2  row  array  is  34  ° .  An  ideal  dipole  is  used  as  the  array  element. 
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Half  power  beamwidth  of  1  row  and  2  rows  array  with  8  elements 


Scan  angle,  degrees 

Figure  15.  Elevation  plane  HPBW  comparison. 

In  order  to  have  the  largest  values  of  gain  in  the  coverage  region,  the  array  can  be 
tilted  back  such  that  its  HPBW  point  touches  the  horizon,  as  shown  in  Figure  16.  The 
range  of  elevation  angles,  as  the  UAV  flies  at  an  altitude  of  300  m,  is  shown  in  Figure  17. 


Figure  16.  Tilting  of  array  for  maximum  elevation  coverage. 


Figure  17.  Elevation  angles  at  various  ranges. 
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The  UAY  is  within  the  HPBW  of  a  one  row  array  for  all  the  angles  shown  in 
Figure  17.  For  the  2  row  array,  the  UAY  is  in  the  HPBW  at  ranges  greater  than  560  m. 
Therefore,  elevation  tracking  would  have  to  be  done  at  close  ranges.  One  advantage  of 
the  2  row  array  is  that  a  null  is  pointed  at  the  ground.  The  null  and  the  low  sidelobe  help 
to  suppress  multipath. 

The  same  Matlab  code  in  Appendix  B  is  used  to  generate  the  azimuth  pattern  of 
the  array  antenna  for  three  different  conditions,  which  are  shown  in  Figure  18.  The  three 
conditions  are  as  follows: 

1 .  1  row  array  at  6S  =90  °  and  <f>s  =90  ° , 

2.  2  row  array  at  6S  =90  °  and  </>s  =90  ° ,  and 

3.  1  row  array  at  0S  =40°  and  (f>s  =90°. 

It  is  noticed  that  the  azimuth  pattern  is  the  same  for  the  first  and  second 

configurations. 
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Figure  18.  Azimuth  (H-plane)  coverage  of  different  arrays. 
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The  elevation  pattern  in  Figure  19  is  generated  for  the  three  different 
configurations  using  the  same  Matlab  codes  in  Appendix  B. 


Figure  19.  Elevation  (E-plane)  coverage  of  different  arrays. 


This  tilting  of  the  array  antenna  will  reduce  the  multipath  effects  and  the 
reflections  from  the  earth.  The  null  region  and  first  sidelobe  of  the  elevation  pattern  for 
the  2  row  array  is  pointed  at  the  ground,  which  reduces  reflections.  Due  to  the  constraint 
on  the  size  of  the  dipole  cards,  the  closest  x-direction  element  spacing  is  10  cm.  Although 
by  adding  8  more  elements,  the  signal-to-noise  ratio  (SNR)  can  be  increased  by  a  factor 
of  8,  this  is  not  an  overwhelming  advantage  as  the  transmitter  power  of  the  UAV  is 
sufficiently  high. 

B.  SUM  AND  DIFFERENCE  BEAM  MATLAB  RESULTS  COMPARED 
WITH  CST  MICROWAVE  STUDIO  RESULTS 

The  results,  using  the  formulas  in  Matlab  codes  for  the  sum  and  difference  beams 
will  be  compared  with  those  generated  by  CST  Microwave  Studio  from  the  thesis  by 
Gezer  [2],  The  reason  for  comparing  the  Matlab  results  with  the  CST  Microwave  Studio 
is  to  verify  that  the  Matlab  results  are  close  to  the  CST  Microwave  Studio  results,  as  CST 
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Microwave  Studio  includes  the  effects  of  the  ground  plane  edges  and  mutual  coupling  on 
the  patterns.  It  also  has  the  exact  dipole  circuit  modeled,  which  gives  the  proper  element 
factor. 


CST  Microwave  Studio  is  a  software  package  for  electromagnetic  analysis  and 
design  in  the  high  frequency  range.  The  Finite  Integration  Technique  (FIT),  a  time- 
domain  numerical  method  based  on  the  integral  form  of  Maxwell’s  equations,  is  used  in 
CST  Microwave  Studio.  The  parameters  of  the  device  that  is  modeled  can  be  defined  in 
detail  in  three  dimensions.  The  program  is  especially  suited  to  the  fast  efficient  analysis 
and  design  of  components  like  antennas,  filters,  transmission  lines,  couplers,  connectors, 
printed  circuit  boards,  resonators,  etc. 


The  simulation  model  used  to  run  the  sum  and  difference  patterns  of  the  array  of 
dipoles  with  scanning  in  the  H-plane  of  the  array  is  shown  in  Figure  20. 


Figure  20.  Top  view  of  the  array  of  dipoles  with  ground  plane  modeled  in  CST 

Microwave  Studio.  From  [2], 
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The  Matlab  codes  in  Appendix  B  are  used  to  generate  the  sum  and  difference 
beam  for  comparison  with  the  CST  Microwave  Studio.  In  order  to  compute  the  difference 
beam,  AFZ  in  Equation  (18)  is  modified  slightly.  Assuming  Nz  is  even,  for  the 
difference  beam: 

Nz/2  Nz 

^fi^cnO-casO,)  (29) 

m- 1  m=Nz/2+\ 

It  can  be  seen  that  the  array  output  of  one-half  of  the  array  is  subtracted  from  the  output 
of  the  other  half. 

The  H-plane  comparison  results  of  the  difference  and  sum  beam  for  a  scan  of  0  ° 
from  broadside  are  shown  in  Figures  21  and  22,  respectively.  The  two  methods  are  in 
good  agreement,  except  at  very  wide  angles. 


AFZ  = 


z 


JPzm  (cos  cos  9S ) . 


z 


Figure  2 1 .  Azimuth  (H-plane)  difference  beam  comparison  at  scan  of  0  °  from 

broadside. 
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Figure  22.  Azimuth  (H-plane)  sum  beam  comparison  at  scan  of  0  °  from  broadside. 

The  H-plane  comparison  results  of  the  difference  and  sum  beam  for  a  scan  of  40  ° 
from  broadside  are  shown  in  Figures  23  and  24,  respectively.  Again,  the  two  methods  are 
in  good  agreement. 


Figure  23.  Azimuth  (H-plane)  difference  beam  comparison  at  scan  of  40  °  from 

broadside. 
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Figure  24.  Azimuth  (H-plane)  sum  beam  comparison  at  scan  of  40  °  from  broadside. 

The  E-plane  comparison  of  results  is  shown  in  Figure  25.  For  the  Matlab 
calculation,  both  a  half-wave  dipole  pattern  and  ideal  dipole  pattern  are  shown,  as 
described  in  Appendix  B. 


Figure  25.  Elevation  (E-plane)  pattern  comparison. 
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c. 


SUMMARY 


The  one  row  array  is  chosen  primarily  because  the  UAV  will  be  within  the  HPBW 
for  all  ranges.  Therefore  elevation  tracking  is  not  necessary.  Although  the  SNR  can  be 
increased  by  8  times  by  adding  a  row  of  8  more  elements,  this  will  add  to  the  cost  of 
equipment  and  the  UAV  transmitter  power  is  sufficiently  high.  The  azimuth  and 
elevation  patterns  of  the  array  were  examined.  The  sum  and  difference  beams  for  the 
Matlab  simulations  are  in  close  agreement  with  the  CST  Microwave  Studio  simulations. 
The  printed  circuit  dipole  element  was  found  to  have  an  E-plane  pattern  HPBW  between 
that  of  an  ideal  dipole  and  a  half-wave  dipole. 
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IV.  FREQUENCY  MODULATION  (FM)  DEMODULATOR 


A.  INTRODUCTION 

Apart  from  tracking  the  UAV,  there  is  a  requirement  to  recover  the  information 
from  the  signals  transmitted  by  the  UAV.  The  signal  of  interest  will  be  either  a  NTSC 
video  signal  or  audio  signal,  that  is  frequency  modulated  (FM)  onto  a  2.4  GHz  carrier. 
The  modulation  characteristics  of  the  FM  video  and  audio  signals  are  shown  in  Table  1 . 


Parameters 

Bandwi  dth 

FM 

Modulation 

Modulation 

Index 

Video  and 
Audio 

6  MHz 

Wideband 

>1.0 

Audio 

22  kHz 

Narrowband 

<1.0 

Table  1.  Parameters  of  FM  video  and  audio  signals.  From  [22], 


UAV  surveillance  has  attempted  to  take  advantage  of  cheap  wireless  video 
transmitters  or  receivers  that  operate  at  2.4  GHz  in  the  ISM  band  (license  free  band). 
Typical  low  power  transmitters  have  power  ratings  in  the  range  of  10-100  mW,  although 
more  powerful  devices  are  available.  For  UAV  applications,  longer  ranges  are  possible 
using  higher  power.  The  2.4  GHz  wireless  video  transmitters,  which  can  transmit  a 
broadcast  quality  picture,  generally  use  a  FM  modulation  method  for  audio/video 
transmission.  One  thing  to  note  is  that  2.4  GHz  wireless  video  systems  are  not  the  only 
equipment  that  uses  this  frequency  band.  For  example  WLAN  systems  and  microwave 
ovens  operate  at  the  same  frequency  band  and  can  cause  interference  in  the  operation  of  a 
2.4  GHz  video  link. 

The  output  of  each  array  element  is  an  FM  signal  with  a  2.4  GHz  carrier.  These 
signals  need  to  be  weighted  and  summed  for  beamforming,  and  then  demodulated  to  get 
the  video  or  audio.  The  beamforming  can  be  done  at  the  carrier  frequency  or  at  baseband. 
Either  way,  it  is  necessary  that  the  incident  wave  phase  change  between  elements  is 
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retained  for  beamforming.  There  are  several  approaches  to  demodulating  the  FM  signal  at 
the  antenna.  Table  2  shows  the  types  of  FM  detectors  that  are  used  for  the  demodulation 
of  FM  signals. 


Name  of 
detector 

Type  of 
detector 

Circuit  used 

Frequency 

discriminator 

Non- coherent 

Quadrature  demodulator  (Bilotti) 
Foster-Seeley  detector 

Balanced  frequency  discriminator 

Phase-locked 
loop  (PLL) 

Coherent 

Analog  phased-locked  loop 

Digital 

sampling 

Coherent 

Analog-to-Digital  Converter 

(ADC) 

Table  2.  Types  of  FM  detectors.  From  [24] 


The  ADC  for  demodulation  of  the  FM  signals  (also  known  as  the  digital 
demodulator),  at  some  point,  converts  the  analog  signal  into  a  digital  representation  prior 
to  demodulation.  The  most  direct  approach  is  to  sample  the  received  signal  out  of  the 
antenna,  but  this  requires  a  very  high-speed  ADC.  If  ten  times  oversampling  is  used 
(which  is  common),  the  sampling  frequency  would  have  to  be  10x2.4  GHz  =  24  GHz. 
The  fastest  ADCs  today  are  in  the  2-4  GHz  range.  Two  common  approaches  to 
implement  this  are  to  sample  at  a  lower  intermediate  frequency  (IF),  e.g.,  around  30  MHz, 
or  to  mix  the  signal  down  to  baseband,  where  I/Q  components  are  required  to  retain  all 
the  information  in  the  signal.  The  Foster-Seely  detector  and  balanced  frequency 
discriminator  convert  FM  into  amplitude  modulation  (AM)  and  use  an  AM  detector  to 
recover  the  message  signal  from  the  converted  waveform.  The  quadrature  demodulator 
converts  FM  into  both  phase  modulation  (PM)  and  AM  and  either  type  of  detector  can  be 
used  to  recover  the  message  signal.  This  is  referred  to  as  Bilotti’s  method. 

A  local  oscillator  (LO)  signal  is  necessary  for  downconverting  from  2.4  GHz  to 
baseband.  The  LO  signal  can  be  provided  by  an  oscillator  at  the  antenna,  or  by  extracting 
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(recovering)  the  carrier  from  the  received  signal.  Each  of  these  methods  has  advantages 
and  disadvantages.  A  phase-locked  loop  (PLL)  is  one  way  to  recover  the  carrier.  A  PLL 
is  an  electronic  control  system  that  generates  a  signal  that  has  a  fixed  relation  to  the  phase 
of  a  "reference"  signal.  A  phase-locked  loop  circuit  responds  to  both  the  frequency  and 
the  phase  of  the  input  signals,  automatically  raising  or  lowering  the  frequency  of  a 
controlled  oscillator  until  it  is  matched  to  the  reference  in  both  frequency  and  phase  [25]. 
If  the  analog  signal  processing,  prior  to  sampling  into  the  digital  domain,  is  done  with 
phase-locked  sources,  the  digital  demodulator  can  support  either  coherent  or  non¬ 
coherent  demodulation.  This  is  a  key  potential  benefit  of  the  digital  solution  [26], 

Several  of  the  methods  for  the  demodulation  of  the  FM  signals  will  be  examined 
to  assess  their  applicability  to  the  tracking  array. 

B.  LOCAL  OSCILLATOR  AND  DEMODULATION  APPROACHES 


Three  FM  demodulation  methods  are  briefly  explained  in  this  section.  A  closely 
related  issue  is  how  the  local  oscillator  signal  is  generated.  The  LO  signal  can  be 
extracted  from  the  incoming  signal  or  generated  locally.  Both  of  these  methods  may 
require  feedback  to  keep  the  LO  locked  on  the  incoming  carrier. 


1.  Phase  Locked  Loop 


A  PLL  is  a  control  circuit  that  is  capable  of  acquiring  a  signal  and  locking  on  to  it. 
The  PLL  system  consists  of  three  parts:  phase  comparator,  low-pass  filter  (LPF),  and 
voltage-controlled  oscillator  (YCO).  All  parts  are  connected  to  form  a  closed-loop 
frequency-feedback  system,  as  shown  in  Figure  26. 


Input 

Signal 


Voltage 


Baseband 

Signal 


Figure  26.  PLL  Block  Diagram.  From  [27], 
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With  no  signal  input  applied  to  the  PLL  system,  the  error  voltage  at  the  output  of 
the  phase  comparator  is  zero.  The  voltage,  Vc!  (t) ,  from  the  LPF  also  is  zero,  which  causes 

the  VCO  to  operate  at  a  set  frequency,  f0,  called  the  center  frequency.  When  an  input 

signal  is  applied  to  the  PLL,  the  phase  comparator  compares  the  phase  and  frequency  of 
the  signal  input  with  the  VCO  frequency.  It  then  generates  an  error  voltage  proportional 
to  the  phase  and  frequency  difference  of  the  input  signal  and  the  VCO.  The  error  voltage, 
Ve(t) ,  is  filtered  and  applied  to  the  control  input  of  the  VCO.  V(/(t)  varies  in  a  direction 

that  reduces  the  frequency  difference  between  the  VCO  and  signal-input  frequency. 
When  the  input  frequency  is  sufficiently  close  to  the  VCO  frequency,  the  closed- loop 
nature  of  the  PLL  forces  the  VCO  to  lock  in  frequency  with  the  signal  input;  i.e.,  when 
the  PLL  is  in  lock,  the  VCO  frequency  is  identical  to  the  signal  input,  except  for  a  finite 
phase  difference.  The  range  of  frequencies  over  which  the  PLL  can  maintain  this  locked 
condition  is  defined  as  the  lock  range  of  the  system.  The  lock  range  always  is  larger  than 
the  band  of  frequencies  over  which  the  PLL  can  acquire  a  locked  condition  with  the 
signal  input.  This  latter  band  of  frequencies  is  defined  as  the  capture  range  of  the  PLL 
system  [27]. 

2.  Local  Oscillator  With  No  Feedback 


Using  an  independent  LO  with  no  feedback  is  similar  to  a  PLL  in  the  free  running 
mode,  as  shown  in  Figure  27.  Ideally  the  VCO  frequency  should  be  exactly  the  same  as 
the  carrier  frequency,  fc .  If  not,  then  the  reference  for  frequency  demodulation  will  be 


off.  If  the  carrier  frequency  of  the  transmitter  is  stable  and  accurate  and  the  VCO 
frequency  is  also  stable  and  accurate,  then  the  baseband  signal  can  be  demodulated.  In  the 


next  two  sections,  Matlab  simulations  are  demonstrated  to  estimate  how  much  frequency 


error  between  the  input  signal  and  VCO  can  be  tolerated. 


Input 

Signal 


Baseband 

Signal 


Figure  27.  Stand-alone  VCO  for  a  local  oscillator. 
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3. 


Bilotti’s  Quadrature  Demodulator3 


The  modulation  of  an  FM  signal  is  the  deviation  of  a  carrier  from  its  nominal 
frequency.  To  demodulate  this  signal,  its  frequency  deviation  is  converted  to  either  phase 
or  amplitude  and  the  change  in  phase  or  amplitude  is  detected.  In  the  quadrature 
demodulator  (also  called  Bilotti’s  method),  the  modulated  carrier  is  passed  through  an  LC 
tank  circuit  that  shifts  the  signal  by  90  0  at  the  center  frequency.  The  direction  of 
frequency  deviation  will  cause  the  phase  shift  to  be  either  greater  or  smaller  than  90  ° . 
The  LC  tank  circuit  simultaneously  converts  the  FM  to  amplitude  modulation  (AM),  so 
an  amplitude  detector  could  also  be  used. 

Bilotti's  quadrature  demodulator  uses  a  phase  shift  network  and  a  phase  detector, 
as  shown  in  Figure  28  [29]. 


Figure  28.  Quadrature  demodulator  block  diagram.  From  [29]. 

The  signal  out  of  the  antenna  is  vin .  The  phase  detector  compares  the  phase  of  the 
input  signal  ( v\ )  to  v2 ,  the  signal  generated  by  passing  v\  through  a  phase  shift  network. 
This  phase  shift  network  includes  an  LC  tank  ( L,Rp  and  C'p )  and  a  series  reactance 

( Cs ).  The  network  gives  a  frequency-sensitive  phase  shift.  The  phase  detector  discussed 
is  the  bipolar,  double-balanced  multiplier  popularized  by  Bilotti.  The  output  of  the 
multiplier  ( i0 )  is  fdtered,  which  results  in  a  DC  level  that  changes  as  the  input  frequency 
changes. 


3  Much  of  this  material  in  section  A  and  B  is  from  [29], 
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To  derive  the  transfer  function  of  the  quadrature  demodulator,  the  phase  shift 
network  is  drawn  as  a  small-signal  circuit  model  in  Figure  29. 


Figure  29.  Small  signal  model  of  the  quadrature  phase-shift  network.  From  [29], 


The  impedance  (Zp)  of  the  parallel  combination  of  L ,  Rp  and  Cp  is: 


/  \  RvLs 

Zp(s)  = - 1 -  (30) 

CpRpLs1  +Ls  +  Rp 

where  5  =  jco.  The  ratio  of  v2  over  vj  is  the  ratio  of  impedances  Zp  (.v)  over 
(Zp  (5)  + 1  /  sCs ) .  Simplifying  this  ratio  gives: 


V2  = _ CsLs2 _ 

Vl  (Cp  +Cs)Ls2  +-^—  s  +  1 


(31) 


Equation  (31)  is  used  in  the  Matlab  code  in  Appendix  C  to  generate  the  plot  of  the  phase 
angle  versus  the  frequency.  The  resonant  frequency  con  of  this  filter  is: 


a"  Jl(Cs+Cp) 


(32) 


40 


C.  SIMULATION  RESULTS 

In  this  section,  simulation  results  are  presented  for  several  FM  demodulation 
techniques. 

1.  Quadrature  Detector  (Bilotti’s  Method) 

In  Bilotti’s  quadrature  detector,  the  output  can  be  applied  to  either  an  AM  or  PM 
detector.  In  the  present  application,  the  mixer  and  filter  in  Figure  28  are  incorporated  in 
the  AD8347.  The  output  iQ  is  in  I/Q  format.  The  demodulation  is  done  in  a  computer 

using  the  sampled  I/Q  data.  Therefore,  either  AM  or  PM  could  be  used  to  recover  the 
video  signal. 

Equation  (31)  and  the  Digikey  catalog  [30]  are  used  to  compute  and  select  the 
appropriate  values  for  the  capacitances  Cp  and  Cs  and  inductance  L  to  obtain  a  carrier 

frequency  of  2.4  GHz.  Using  off  the  shelf  values  of  Cp  and  Cs ,  both  are  5  pF.  From 

available  off  the  shelf  values  of  L ,  45  nH  gives  the  exact  carrier  frequency  of  2372.542 
MHz.  These  lumped  element  values  are  used  in  the  simulations. 

Using  the  Matlab  codes  in  Appendix  C  and  varying  the  values  of  Rp ,  various 

plots  for  the  6  MHz  video  bandwidth  are  obtained  as  shown  from  Figures  30  to  33. 
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Frequency,  MHz 


Figure  30.  Magnitude  of  —  plot  at  Rv  =560  Q ,  1000  Q  and  1500  Q . 

V1 


Frequency,  MHz 


Figure  31.  Angle  of  —  plot  at  Rp  =560  Q ,  1000  Q  and  1500  Q . 
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Frequency,  MHz 


Figure  32.  Magnitude  of  —  plot  at  Rv  =2700  Q ,  3900  Q ,  5600  Q  and  10  kQ . 

vl 
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It  can  be  observed  that  as  Rp  increases,  more  phase  angles  can  be  covered  within 

the  6  MHz  bandwidth  about  the  carrier  frequency  of  2372  MHz.  However,  the  plots 
become  more  and  more  non-linear.  The  characteristics  are  summarized  in  Table  3  for  the 
6  MHz  bandwidth  at  about  2372  MHz. 


Angle  range,0 

Behaviour 

560 

80  to  104 

Linear 

1000 

71  to  114 

Linear 

1500 

64  to  123 

Linear 

2700 

51  to  140 

Non-linear 

3900 

40  to  150 

Non-linear 

5600 

29  to  157 

Non-linear 

10000 

20  to  166 

Non-linear 

Table  3.  Summary  of  the  plots  in  Figures  30  to  33. 


If  the  computer  processor  uses  AM,  then  a  large  dynamic  range  in  amplitude  is 
best.  However,  if  PM  is  used  then  a  linear  phase  curve  over  a  large  range  of  angles  is  best. 
Using  the  Matlab  codes  in  Appendix  C  and  varying  the  values  of  Rp ,  various 

plots  for  the  22  kHz  audio  bandwidth  are  obtained,  as  shown  in  Figures  34  to  37. 


x  104 


Figure  34.  Magnitude  of  —  plot  at  Rp  =510  kQ ,  604  kQ ,  715  kQ  and  806  kQ . 

V1 
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Figure  35.  Angle  of  —  plot  at  RD  =510  kQ ,  604  kQ ,  715  kQ  and  806  kQ . 

vl 


x  10s 


Figure  36.  Magnitude  of  —  plot  at  R„=  1  MQ ,  2  MQ ,  3  MQ  and  4.02  MQ  . 

V1 
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Figure  37.  Angle  of  —  plot  at  Rp=  1  MQ,2MQ,3MQ  and  4.02  M  Q . 

It  is  noted  that  as  Rp  increases  from  510  kQ  to  4.04  MQ,  the  higher  values 

cover  a  larger  range  of  phases,  but  the  non-linear  characteristic  becomes  more  obvious. 
For  the  22  kHz  bandwidth,  Rp  is  chosen  to  be  1  MQ . 

The  quadrature  circuit  could  be  implemented  in  the  digital  array,  as  shown  in 
Figure  38.  The  reference  signal  can  be  obtained  using  a  special  auxiliary  antenna,  or  it 
could  be  tapped  off  of  one  of  the  elements.  However  it  is  done,  the  same  LO  source  must 
be  used  for  all  elements. 
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Figure  38.  Digital  array  using  a  quadrature  detector. 


Figure  39  shows  a  Matlab  Simulink  model  for  the  array  in  Figure  38,  up  to  the 
computer  block.  A  two-tone  signal,  one  sine  wave  at  900  Hz  and  another  at  1700  Hz,  is 
used  as  a  source  that  is  converted  to  FM.  The  simulation  operates  at  50  kHz,  rather  than 
2.4  GHz,  to  greatly  reduce  the  simulation  time.  The  frequency  deviation  of  the  FM  is  set 
at  1000  Hz,  which  gives  a  0.5  modulation  index,  assuming  a  maximum  modulating 
frequency  of  an  FM  signal  is  2000  Hz.  A  sampling  frequency,  which  is  10  times  the 
simulation  frequency,  is  used  to  obtain  a  more  refined  simulation  result.  The  same 
sampling  frequency  will  be  used  for  subsequent  models  in  this  section.  The  RLC  circuit 


is  modeled  by  its  transfer  function,  and  its  resonant  frequency  was  set  to  50  kHz. 


Scope 


Figure  39.  Simulink  model  of  Bilotti’s  quadrature  demodulator. 
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The  FFT  of  the  source  and  recovered  signals  are  shown  in  Figure  40.  The 
frequency  components  of  the  two-tone  signal  are  clearly  visible  in  the  recovered  signal. 
This  shows  that  the  source  signal  can  be  recovered  using  Bilotti’s  method. 


Figure  40.  FFT  of  the  recovered  and  source  signals  in  Bilotti’s  simulink  model. 

2.  Local  Oscillator  With  No  Feedback 


Figure  41  shows  the  block  diagram  for  an  array  with  an  independent  LO.  The  LO 
will  synchronize  the  AD8347  demodulators.  As  this  is  a  method  without  feedback, 
simulations  are  done  to  consider  the  effects  of  a  phase  drift  in  the  LO.  The  raw  I  and  Q 
components  from  the  AD8347  demodulators  are  processed  to  perform  the  beamforming 
and  demodulate  the  FM  video  signal.  The  carrier  frequency  of  50  kHz  with  a  1  kHz 
deviation  is  used  in  the  simulation  for  the  modulator  for  transmitting  the  two-tone  signals. 
Figure  42  is  the  Simulink  model  of  Figure  41. 
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Figure  41.  Open  loop  digital  quadrature  detector. 


Spectrum 

Scopel 

Figure  42.  Simulink  model  of  open  loop  digital  quadrature  detector. 

Figure  43  shows  the  input  and  demodulated  signals  in  the  time  domain.  Figure  44 
shows  their  spectra  in  the  frequency  domain.  The  recovered  signal  matches  the  source 
signal  in  both  the  time  scope  and  the  FFT  spectrum  scope. 
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Two-tone  signal  Magnitude,  volts 


Time,  seconds 

Figure  43.  Time  scope  of  the  open  loop  digital  quadrature  detector  in  Simulink. 


Figure  44.  FFT  spectra  of  the  source  and  recovered  signals. 
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3.  Digital  PLL  Quadrature  Detector 


Figure  45  shows  how  the  PLL  is  integrated  into  the  digital  array  to  perform  the 
beamforming  and  the  FM  demodulation.  The  output  I  and  Q  components  are  converted  to 
the  digital  format  to  be  processed  in  the  computer.  In  the  computer,  the  digital  I  and  Q 
signals  are  low  pass  filtered  before  performing  the  beamforming  and  FM  demodulation 
The  error  signal  is  computed  and  fed  to  the  variable  control  oscillator  (YCO)  to  adjust  the 
VCO  frequency.  The  same  PLL  is  used  for  mixing  all  the  elements  from  the  other 
channels. 


FM  signal 


Video 


Figure  45.  Digital  PLL  quadrature  detector. 


The  Simulink  model  in  Figure  46  was  used  to  simulate  the  digital  PLL  quadrature 
detector  in  Figure  45.  The  carrier  frequency  of  50  kHz  with  a  1  kHz  deviation  is  used  for 
the  modulator  for  transmitting  the  two-tone  signals.  The  VCO  with  the  following 
parameters  is  inserted  in  the  feedback  loop:  output  amplitude  of  3  V,  quiescent  frequency 
of  50  kHz,  input  sensitivity  of  500  Hz/V  and  a  sampling  frequency  of  500  kHz. 
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Figure  46.  Simulink  model  of  digital  PLL  quadrature  detector. 


Figure  47  shows  the  input  and  demodulated  signals  in  the  time  domain.  Figure  48 
shows  their  spectra  in  the  frequency  domain.  It  is  noticed  that  the  source  signal  matches 
the  recovered  signal  for  the  digital  PLL  quadrature  detector. 
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Figure  47.  Time  scope  of  the  digital  PLL  quadrature  detector  in  Simulink. 


52 


20 


Figure  48.  FFT  spectrum  of  the  digital  PLL  quadrature  detector  in  Simulink. 


D.  PHASE  ERRORS  AND  LO  FREQUENCY  DRIFT 

There  is  a  need  to  study  the  effects  of  the  frequency  drift  in  the  LO  and  the 
relative  phase  errors  in  the  array  channels.  The  relative  phase  errors  cause  beamforming 
errors,  which  lead  to  tracking  errors.  In  this  section,  the  effects  of  phase  errors  of  a  4- 
channel  array  modeled  in  Simulink  are  shown.  The  frequency  drift  in  the  LO  is  also 
simulated. 
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1. 


Relative  Phase  Errors  in  the  Array  Channels 


The  Simulink  block  diagram  designed  to  model  the  phase  errors  in  the  array 
channels  is  shown  in  Figure  49  below.  The  carrier  frequency  of  50  kHz  with  a  1  kHz 
deviation  is  used  in  the  simulation  for  the  modulator  for  transmitting  the  two-tone  signals. 
Phase  shifters  are  introduced  to  the  open  loop  digital  quadrature  detector  Simulink  model. 
Different  delay  line  sizes  are  introduced  to  the  other  three  channels  to  model  the  effects 
of  phase  shifters  on  the  incoming  signal.  Ideally,  the  delays  cause  the  sum  output  to  have 
a  maximum  at  the  signal  angle  of  incidence.  The  Simulink  results,  shown  in  Figure  50, 


indicate  that  the  signal  could  be  recovered  when  phase  shifter  errors  are  introduced. 


Figure  49.  Simulink  model  for  phase  errors  simulation  analysis. 
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Figure  50.  Simulation  results  for  phase  errors  in  a  4-channel  array. 

2.  LO  Frequency  Drift 

The  frequency  drift  in  the  LO  is  simulated  in  the  Simulink  model,  shown  in 
Figure  51,  at  a  carrier  frequency  of  5  MHz  to  analyze  the  effects  to  the  recovered  signal. 
A  frequency  offset  of  100  Hz  is  introduced  to  the  transmitter  frequency. 
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Spectrum 


Figure  5 1 .  Simulink  model  for  LO  frequency  drift  analysis. 


The  simulation  result  is  shown  in  Figure  52. 


Figure  52.  Simulation  result  at  5  MHz  with  100  Hz  frequency  drift. 
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E.  SUMMARY 


Simulation  tools  are  developed  for  array  analysis  of  the  three  FM  demodulation 
techniques:  Bilotti’s  method,  local  LO  with  no  feedback  and  digital  phase  locked  loop 
quadrature  detector. 

The  LO  with  no  feedback  is  preferred  as  it  uses  the  least  circuitry  with 
satisfactory  performance  for  the  recovery  of  the  two-tone  signal.  Beamforming  can  also 
be  performed  easily,  as  the  raw  I  and  Q  data  is  preserved.  This  technique  also  eases  the 
computer  processing.  Relative  phase  errors  in  the  array  channels  and  LO  frequency  drift 
are  introduced  into  the  Simulink  model.  The  simulation  results  indicate  that  the  recovery 
of  the  two-tone  signal  can  still  be  achieved.  Hence,  the  LO  with  no  feedback  technique 
can  be  used. 
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V.  SYSTEM  ARCHITECTURE 


This  chapter  discusses  the  hardware  necessary  to  build  the  tracking  array.  To  keep 
the  cost  down,  commercial-off-the-shelf  hardware  is  used. 

A.  DESCRIPTION  OF  TRACKING  ARRAY  HARDWARE 


A  general  block  diagram  of  a  digital  receiver  channel  is  shown  in  Figure  53.  The 
received  FM  video  or  audio  signal  undergoes  the  following  processes: 

1 .  Amplified  by  the  low  noise  amplifier  (LNA)  and  then  fed  to  the  RF  input 
of  an  AD 8 3 47  quadrature  demodulator. 

2.  Downconverted  to  baseband  by  the  AD8347  quadrature  demodulator, 
operating  from  800  to  2700  MHz  with  a  modulation  bandwidth  of  about  50 
MHz.  The  output  signals  are  in  I/Q  format. 

3.  Sampled,  filtered  and  then  sent  to  the  processor. 

4.  Processing  of  the  detected  signals  is  done  to  form  the  beam  by  combining 
the  I/Q  signals  from  all  receive  channels.  The  FM  is  removed  and  the  video 
signal  is  recovered. 

5.  A  common  local  oscillator  is  used  to  synchronize  the  devices  to  avoid 
phase  and  drifting  errors  between  the  elements. 


These  processes  are  illustrated  in  Figure  53. 


Video 

Signal 


Figure  53.  Array  architecture  for  one  channel. 
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The  dipoles,  LNA,  demodulator,  local  oscillator,  power  supplies,  power 
adjustment  circuit  and  signal  processing  are  discussed  briefly,  and  follow  the  description 
in  Gezer  [2], 

1.  Dipoles 

The  dipoles  were  designed  in  the  CST  Microwave  Studio  by  Professor  David  C. 
Jenn  and  fabricated  by  Cirexx  Corporation.  They  operate  at  2.4  GHz.  The  return  loss  of 
the  dipoles  was  measured  and  found  to  be  at  least  15  dB.  Reference  [8]  discusses  the 
dipole  element  return  losses  more  in  detail. 

2.  Low  Noise  Amplifier  (LNA) 

The  LNA  amplifies  the  weak  signals  captured  by  the  dipoles.  The  LNA  has  a  low 
noise  figure  (2.0  dB)  and  high  gain  (24  dB),  which  improves  the  performance  of  the  array. 
Because  the  overall  noise  figure  is  dominated  by  the  first  few  stages,  the  total  noise  of  the 
cascaded  stages  is  reduced  by  the  LNA.  Hence,  the  fact  that  the  LNA  has  a  very  low 
noise  figure  is  a  key  factor  in  signal  detection. 

The  LNA  must  be  located  very  close  to  the  dipole  in  order  to  minimize  the  losses 
in  the  feed  line.  For  our  setup,  the  LNAs  are  connected  to  the  dipoles  with  a  barrel  to 
minimize  the  cable  loss.  This  SNR  will  be  verified  later  in  Section  B.  The  LNA-2700  by 
RF  Bay  Inc.  (shown  in  Figure  54)  was  selected  as  the  LNA  because  of  its  reasonable 
price  as  well  as  its  low  noise  figure. 


Figure  54.  LNA  2700  by  RF  Bay  Inc.  From  [2], 
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3. 


Demodulator 


The  demodulator  downconverts  the  signal  from  the  carrier  frequency  to  the 
baseband.  The  demodulator  selected  for  the  array  is  the  AD8347  Quadrature 
Demodulator  manufactured  by  Analog  Devices.  Figure  5  shows  one  of  these 
demodulators  used  in  the  system.  This  COTS  demodulator  operates  within  the  frequency 
range  0.8  through  2.7  GHz  and  requires  supply  voltages  between  2.7  and  5.5  V.  The 
voltage  YGIN  determines  the  amplification  of  the  received  signal.  Reference  [31] 
includes  a  very  detailed  spec  sheet  of  the  AD8347  EVAL  Quadrature  Demodulator. 

4.  Local  Oscillator 

The  LO  is  designed  and  manufactured  by  Z-Communications.  The  LO  signal  at 
2.4  GHz  is  generated  by  a  Y800ME10  oscillator  on  the  board  and  fed  into  the 
demodulators. 


5.  Power  Supplies  and  Power  Adjustment  Circuit 

The  two  linear  power  supplies  used  for  the  array  are  model  HP  6255  A  by 
Hewlett  Packard.  The  power  supplies  provide  +5  and  +12  V.  The  +12  V  is  applied  to 
VGIN  of  the  demodulator  after  being  adjusted  with  a  VGIN  voltage  regulator,  designed 
by  Bob  Broadston  from  NPS.  The  LO  and  LNAs  are  provided  with  +12  V  by  the  power 
supply,  which  has  maximum  currents  of  6.00  A,  2.50  A  for  the  voltages  of  5  V  and  12  V, 
respectively. 

6.  Analog-to-digital  Converters 

The  ADCs  sample  the  baseband  I/Q  signals  out  of  the  AD8347.  The  minimum 
sampling  rate  is  set  by  the  baseband  bandwidth,  which  is  approximately  6  MHz  for  the 
FM  video  signal.  For  5  times  oversampling,  this  would  be  60  MHz.  For  the  audio  signal, 
it  would  be  about  110  kHz. 

Candidate  hardware  for  the  audio  case  is  the  National  Instruments  Compact 
Reconfigurable  I/O  (cRIO)  product  line.  The  reconfigurable  chassis  is  the  heart  of  the  NI 
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cRIO  system  and  contains  a  Field  Programmable  Gate  Array  (FPGA).  The  cRIO  core  is 
powered  by  a  user-programmable  FPGA  configured  by  a  special  software  called 
LabVIEW,  that  serves  as  a  high-speed  interface  to  the  I/O  modules.  The  LabVIEW 
FPGA  Module  is  used  to  program  the  FPGA  on  the  cRIO  reconfigurable  chassis.  FPGAs 
are  silicon  chips  with  unconnected  logic  gates.  The  flexible  design  of  FPGA  enables  the 
user  to  target  any  device  compatible  with  NI  cRIO  and  to  customize  hardware  functions 
for  its  system.  It  is  possible  to  execute  multiple  tasks  simultaneously  using  the  LabVIEW 
FPGA  module. 

NI  provides  different  I/O  modules  for  different  purposes.  The  four  analog  output 
signals  from  the  demodulator  are  input  to  an  analog  input  (A/I)  module  to  be  sampled. 
They  are  then  converted  into  digital  signals  for  signal  processing  and  analysis.  The  A/I 
module  selected  for  the  audio  case  is  the  NI  cRIO-9215.  The  NI  cRIO-9215,  shown  in 
Figure  55,  provides  4  simultaneously  sampled  analog  inputs  with  a  rate  of  100  kS  and  a 
resolution  of  16  bits. 


Figure  55.  NI  cRIO-9215  Analog  Input  Module.  From  [32], 

The  NI  cRIO-9215  analog  input  modules  are  connected  directly  to  the  FPGA 
hardware.  This  FPGA  hardware  is  the  reconfigurable  chassis  NI  cRIO-9104  with  the  real¬ 
time  embedded  controller  cRIO-9004.  The  cRIO-9104  chassis  has  8  slots  and  a  RIO 
FPGA  core  with  a  RAM  of  196  kB.  The  FPGA  core  has  an  individual  connection  to  each 
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NI  cRIO-9215  analog  input  module  and  can  read  information  from  each  module.  The 
module  operation  is  precisely  synchronized  with  25  ns  resolution. 


B.  RECEIVER  MODULE  ANALYSIS 

This  section  presents  an  estimate  of  the  operating  range  based  on  the  minimum 
power  required  by  the  demodulator  (AD8347). 

1.  Sensitivity  and  Dynamic  Range 

The  maximum  RF  power  that  is  allowed  at  the  receiver  input  (without  saturating 
it)  is  limited  by  the  demodulator’s  1  dB  compression  point  (IPldB).  From  [9],  the 
demodulator’s  IPldB  «  -  22  dBm  with  no  Automatic  Gain  Control  (AGC).  With 
reference  to  Figure  53,  the  LNA  gain  is  24  dB  and  the  cable  losses  on  the  two  sides  of  the 
LNA  are  0.2  dB  and  0.4  dB.  The  maximum  RF  power  that  is  allowable  to  the  receiver 
antenna  is: 

P  =  -22  +  0.2  -24  +  0.4=  -45.4  dBm 

Mnax 

P  gives  the  closest  range  without  over  driving  the  demodulator. 

'max  ^  ^ 

The  minimum  RF  power  that  is  required  at  the  receiver  module  for  operation 
depends  on  the  demodulator’s  sensitivity.  From  [9],  the  demodulator  mixer  output 
voltage  is  linear  with  RF  input  power  from  -  65  dBm  to  -  5  dBm.  Hence,  with  reference 
to  Figure  53,  the  minimum  RF  power  at  the  receiver  antenna  is: 

P.  =  -65  +  0.2  -  24  +0.4  =  -  88.4  dBm 

Mnin 

= -118.4  dBW=  1.445  x  1012W  (33) 

Hence,  the  dynamic  range  of  the  receiver  module  is  -45.4-  (-88.4)  =  43  dB. 

2.  Maximum  Range 

Consider  the  distributed  array  block  diagram  in  Figure  56  where 

Rn  =the  distance  between  the  UAV  and  the  nth  receive  element, 
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Pt  =  transmitted  power  from  the  UAV, 

P  =  received  power  at  nth  element, 

Gt  =  gain  of  transmitter  antenna  of  UAV,  and 

G0  =  gain  of  receiver  antenna  of  n‘h  element  (assume  equal  gains,  G0  =  G0 ). 

UAV 


Figure  56.  Distributed  array  of  receiver  antennas. 


At  any  one  time,  N  elements  are  receiving  signals  from  the  UAV,  where  N  is  the 
total  number  of  elements  in  the  array,  which  in  our  case  is  N  =  8.  From  Figure  56,  the 
one-way  Friis  equation  is  used  to  determine  the  maximum  operating  range,  Rm,ix .  It  is 

considered  that  there  is  no  polarization  mismatch,  input  mismatch  or  other  system  losses. 
The  one-way  signal  strength  from  a  transmitter  to  a  receiver  element  is: 


P  = 


(4  4K 


(34) 


To  include  receiver  sensitivity  in  Equation  (34),  P  is  used  to  determine  the 

rmin 

maximum  range: 


I 


(35) 
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It  will  be  assumed  that  Pt  is  100  mW  in  this  case.  The  gain  of  each  of  the  receiver 
elements  is  G0  and  calculated  to  be  6.3  =  8  dBi.  Thus,  substitute  Pt  =  100  mW,  Gt  =  1.6, 
Gn  =  6.3,  and  P  =  -88.4  dBm  into  Equation  (35)  to  get  f?  =  8.31  km. 

U  'min  mdA 

C.  SUMMARY 

A  description  of  tracking  array  hardware  was  the  subject  of  this  chapter.  The 
possible  hardware  components  were  briefly  introduced.  From  the  block  diagram  of  the 
receiver  module,  the  minimum  power  for  the  receiver  module  was  used  to  determine  the 
maximum  operating  range  of  the  array  using  the  one-way  Friis  equation. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY  AND  CONCLUSIONS 

As  a  part  of  an  ongoing  project,  this  research’s  aim  was  to  design  a  digital  phased 
array  for  receiving  and  tracking  a  2.4  GHz  FM  video  signal  from  an  unmanned  air 
vehicle.  Tracking  the  signal  from  the  UAV  will  establish  a  data  link  for  data  transfer  such 
as  video,  image  or  audio  between  the  UAV  and  the  ground  station. 

An  array  tradeoff  study  was  done  to  determine  the  benefits  of  the  1  row  and  2  row 
arrays  with  8  elements  per  row  in  terms  of  HPBWs  of  the  azimuth  and  elevation  patterns. 
The  one  row  array  was  chosen  primarily  due  to  its  wide  elevation  pattern  coverage.  This 
eliminates  the  need  for  elevation  tracking,  as  the  UAV  is  within  the  HPBW  for  all  UAV 
ranges  of  interest.  Adding  a  row  of  8  more  elements  to  increase  the  SNR  will  increase  the 
cost  of  equipment.  The  additional  SNR  is  not  required,  as  the  UAV  transmitter  power  is 
sufficiently  high.  The  sum  and  difference  beams  for  the  Matlab  simulation,  based  on  the 
derived  pattern  function  of  the  one  row  array,  are  in  agreement  with  the  CST  Microwave 
Studio  simulations  from  the  previous  project  by  Gezer  [2],  The  printed  circuit  dipole 
element  was  found  to  have  an  E-plane  pattern  HPBW  (80°)  between  that  of  an  ideal 
dipole  (90° )  and  a  half-wave  dipole  (78  ° ). 

Several  approaches  for  the  demodulation  of  FM  signals  were  explored.  The 
considerations  for  the  choice  of  the  approaches  for  the  FM  demodulation  are  minimum 
circuitry,  the  ability  to  perform  beamforming,  recovery  of  the  transmitted  signal  and  the 
susceptibility  to  phase  errors  and  LO  drift.  Three  FM  demodulation  techniques  are 
considered:  Bilotti’s  quadrature  demodulator,  a  local  oscillator  with  no  feedback  and  a 
digital  phase  locked  loop  quadrature  detector.  All  three  use  digital  beamforming  and 
signal  processing  to  demodulate  the  FM. 

The  Simulink  models  of  the  three  FM  demodulation  techniques  show  satisfactory 
performance  for  the  recovery  of  a  two-tone  signal.  The  simulations  are  done  at  a  carrier 
frequency  of  50  kHz,  instead  of  2.4  GHz,  to  greatly  save  on  the  simulation  time.  A  LO 
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with  no  feedback  is  chosen,  as  it  uses  the  least  circuitry  with  satisfactory  performance  for 
the  recovery  of  the  two-tone  signal.  Beamforming  can  also  be  performed  easily,  as  the 
raw  I  and  Q  data  are  preserved.  This  technique  also  eases  the  computer  processing. 
Simulations  were  done  with  phase  errors  introduced  in  the  channels  and  a  frequency 
offset  in  the  LO  without  feedback.  The  Simulink  results  indicate  that  the  recovery  of  the 
two-tone  signal  can  still  be  achieved  with  the  expected  LO  drift. 

Finally,  the  possible  hardware  components  are  briefly  introduced  in  Chapter  Y. 
From  the  block  diagram  of  the  receiver  module,  the  minimum  power  for  the  receiver 
module  was  used  to  determine  the  maximum  operating  range  of  the  array  using  the  one¬ 
way  Friis  equation. 

B.  RECOMMENDATIONS  FOR  FUTURE  WORK 

Simulations  and  studies  should  be  conducted  to  analyze  the  following  effects  on 
the  array  and  demodulation  of  FM  signals. 

1.  Implementation  of  Signal  Processing 

The  software  necessary  to  do  the  beamforming  and  FM  demodulation  must  be 
acquired  or  developed.  The  beamforming  and  tracking  are  relatively  straightforward,  and 
have  been  demonstrated  by  Gezer  [2],  However,  an  acquisition  method  must  be  devised 
to  acquire  the  signal  when  the  antenna  and  processor  are  first  turned  on  (or  when  the 
transmitter  signal  first  appears).  If  the  PLL  method  is  used,  the  processor  must  compute 
the  VCO  error  signal. 

2.  Hardware  Tradeoff  Studies 

An  appropriate  ADC  must  be  selected  and  interfaced  with  the  AD8347 
demodulator  board,  and  controller  software  must  be  written.  The  hardware  must  be 
selected  to  minimize  the  amount  of  processing  and  control  commands,  so  that  real-time 
operation  is  possible. 
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APPENDIX  A.  DIPOLES  OF  ARBITRARY  ORIENTATION4 


The  dipole  is  centered  at  the  origin  as  shown  in  Figure  Al.  The  dipole  direction 
is  given  by  the  vector: 

I  =  xsin 6a  cos,(/)a  +y$>m6a  sin^a  +  zcos0a 

ua  Va  wa 

The  dipole  direction  cosines  are: 

v  direction  cosine:  ua  =  cos  ax  =  sm6a  cos  (j)a 
y  direction  cosine:  va  =  cos  ay  =  sin  0a  sin  (j)a 
z  direction  cosine:  wa  =  cosaz  =  cos  6a 

The  length  of  the  dipole  is  L,  and  for  a  half-wave  dipole  the  current  distribution  is  given 
by  the  following  function  of  the  path  length  variable  l' : 

7(0  =  Im  cos {ki'\  -L12<t’<LI2 

where  Im  is  the  maximum  value  at  the  feed  point.  For  an  ideal  (Hertzian)  dipole  with 
constant  current: 


1(f)  =  Im ,  -L/2<f<L/2,  L«A 

The  direction  cosines  for  the  far  field  observation  point  are: 

x  direction  cosine:  u  =  sin  6 cos  ^ 
y  direction  cosine:  v  =  sin  6  sin  (/> 
z  direction  cosine:  w  =  cos  8 


4  The  materials  for  Appendix  A  is  from  “Arbitary  dipoles”  unpublished  notes  of  David  C.  Jenn. 
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Figure  Al:  Dipole  of  arbitrary  orientation  centered  at  the  origin 
The  far- field  radiation  integral  can  be  written  as: 


E(0,</>) 


jPvJ_e-jPr 

Anr 


-L/2 


where  it  is  understood  that  the  r  component  is  discarded.  For  the  ideal  dipole  the 
integral  is: 


{  I(£')ejPe'^cU'  =  ImL  s  inc 
-L/2 


'  J3Lr*P 
2  , 


'  , m L 


(L  «  X) 


and  the  resulting  spherical  field  components  are: 

E0  =  E(0,  </>)  •  0  =  ~  'f3,1oImL  e-jprj  #  q 


Anr 


e6  =  m  $)'i=  jP?°ImL  e~iPri 


Akv 
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From  the  transform  tables: 


0  =  x  cos  8  cos  <j)  +  y  cos  8  sin  $  -  z  sin  8 
(f>  =  -xsin^  + j>cos^ 


For  a  Hertzian  dipole  along  the  z  axis:  ( 8a  =  (j)a  =  0° , 7  =  z) 


I  •  8  =  z*  8  =  -sind 

/V  /V  /V 

I  •</>  =  z*t/>  =  0 

and  the  fields  reduce  to  the  familiar  results: 


Ee  =  J/3,1°1»iL  e~JPr  sin  Q 


E^=0 


47rr 


For  a  half-wave  dipole  (L  =  XI 2)  the  integral  is: 


LI  2 

J  Im  cos(/?f 1 ')e 
-Lt  2 


Imej Mr*1)  [  j  ftf  (p  •  i)  cos(flf)  +  fl  sin(/3f)\ 


[fM'(r*i)]  +J37 


L!  2 


-L!  2 


After  applying  some  trigonometric  identities,  the  final  result  for  the  far  electric  field  is: 


E(8,(/>)  = 


jnjJEEL  cos^r./ 

lKr  l-(r./)2 


For  a  half-wave  dipole  along  the  z  axis:  Qa  =  <f>a  =  0° ,7  =  z 

r  •  I  =  cos  8  =  w 
<9  •  /  =  -  sin  # 


F  _  J’loL, e  ,r" 

Ee  2nr 

E(j)=  0 


sin  (9 
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For  a  half-wave  dipole  along  the  x  axis:  6a  =  90°  ,<f>a  =  0° ,  /  =  x 


r  •/  =  sin#cos^  =  u 

/y  /y 

6*1  =  COS&COS0 

/ \  /V 

</>•!  =  -sin^ 


En  = 


JVo^m  e 


~jpr 


2  nr 


cos  sin#  cos 
1-sin2  #cos2  (f) 


JVo^m  e 


-jfir 


2nr 


cos  — 


|ySin#cos^ 


1-sin2  6 cos2  (/) 


cos  6  cos  (j> 


sin  (j> 
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APPENDIX  B.  MATLAB  CODE  FOR  AN  ARRAY  OVER  A 

GROUND  PLANE 


This  Matlab  code  was  used  to  compute  the  array  pattern  for  1  row  and  2  row  arrays  and  the  sum 
and  difference  beams,  using  the  formulas  derived  in  Chapter  II. 


linear  array  of  dipoles  with  a  ground  plane 
%  z  is  array  axis;  ground  plane  is  the  xz  plane 
%  y  is  normal  to  ground  plane 
%  option  to  select  dipoles  along  x  or  z 
%  fixed  to  include  Etheta  and  Ephi 

%  two  dimensional  array  (scans  in  both  x  and  z  dimension) 

%  half-wave  dipoles 
clear 

warning  off 

rad=pi/180; 

f=2.4e9; 

wave=3e8/f; 

beta=2*pi/wave; 

dz=0.065; 

Nz=8;  %no.  of  elements  in  a  row 
dx=0.1; 

Nx=l;  %To  define  the  no.  of  rows 
h=0.25*wave; 

%  select  the  dipole  direction:  V  or  'z' 
dipdir-x';  %  dipoles  x  directed  (parallel  case) 

%dipdir=’z';  %  dipoles  z  directed  (collinear  case) 

%  select  the  dipole  type:  'h -halfwave,  ?s -short  (constant  current) 
diptype-s'; 

%diptype-h'; 

%phi=0  cut 

%  phi=0  vary  theta  (x-z  plane) 
phi=90*rad; 

theta=[0: 1:180] *rad;  w=cos(theta); 

thetas=90;  phis=;%theta  defines  the  yz  plane;  phi  defines  the  xy  plane; 

%  x  is  along  the  direction  of  the  short  dipole; 

%y  is  perpendicular  to  the  ground; 

%z  is  perpendicular  to  the  short  dipole  and  parallel  to  the  ground  plane. 
cts=cos(thetas*rad);  u=sin(theta)*cos(phi);  v=sin(theta)*sin(phi); 
cps=cos(phis*rad); 
alphaz=-beta*  dz  *  ct  s ; 
alphax=-beta*dx*cps; 
psiz=beta*dz*w+alphaz; 
psix=beta*dx*u+alphax; 
if  dipdir=- z'  %  colinear  dipoles  (z  directed) 
if  diptype==V 

EF  theta=2  *  sqrt(  1  -w.  A2) .  *  sin(beta*h*  v) ; 

EFphi=0; 

end 

if  diftype=- h’ 

EFtheta=2*sqrt(l-w.A2).*sin(beta*h*v).*cos(pi*w/2)./sqrt(l-w.A2+le-5); 
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EFphi=0; 

end 

end 

if  dipdir=='x’  %  parallel  case  (x  directed) 
if  diptype=- s' 

EF  theta=2  *  cos(phi) .  * w. *  sin(beta*h*  v) ; 

EFphi=-2  *  sin(phi) .  *  sin(beta*h*  v) ; 
end 

if  diptype=- h' 

EFtheta=2*cos(phi).*w.*sin(beta*h*v).*cos(pi*u/2)./(l-u.A2+le-5); 

EFphi=-2*sin(phi).*sin(beta*h*v).*cos(pi*u/2); 

end 

end 

AF=sin(Nx*psix/2)./sin(psix/2).*sin(Nz*psiz/2)./sin(psiz/2); 

caftheta=abs(AF .  *EF  theta)/2/Nx/Nz; 

cafphi=abs(AF.*EFphi)/2/Nx/Nz; 

cafthdb=20  *  log  1 0(caftheta+ 1  e-5 ) ; 

cafphdb=20*logl0(cafphi+le-5); 

caf=sqrt(cafphi.A2+caftheta.A2); 

cafdb=  1 0  *  log  1 0(caf.  A2) ; 

figure(l),clf 

subplot(121) 

polardb(theta/rad,cafthdb,'-') 

title(['Polar  dB  plot  EYtheta,  \phi  ,num2str(phi/rad),'Ao']) 
subplot(122) 

polardb(theta/rad,cafphdb,'-') 

title(['Polar  dB  plot  E_\phi,  \phi=',num2str(phi/rad),'Ao']) 
hold  off 
%theta=90  cut 

clear  theta  phi  AF  caf  cafthdb  cafphdb  EFth  EFph  psi  u  v  w  cts  ct 
%  theta=90  vary  theta  (x-y  plane) 
phi=[0:l:180]*rad; 
theta=89.99*rad;  ct=cos(theta); 

%  note  to  always  mantain  the  cut  theta  at  89.99  or  90.01 
u=sin(theta)  *  cos(phi) ; 

cts=cos(thetas*rad);  v=sin(theta).*sin(phi);  w=ones(l  ,length(v))*ct; 
alphaz=-beta*  dz  *  ct  s ; 

psiz=(beta*dz*cos(theta)+alphaz)*ones(length(phi),l)'; 
psix=(beta*dx*u+alphax); 
if  dipdir=- z'  %  colinear  dipoles  (z  directed) 
if  diptype=- s' 

EF  th=2  *  sqrt(  1  -w.  A2) .  *  sin(beta*h*  v) ; 

EFph=0; 

end 

if  diftype=='h' 

EFth=2*sqrt(l-w.A2).*sin(beta*h*v).*cos(pi*w/2)./sqrt(l-w.A2+le-5); 

EFph=0; 

end 

end 

if  dipdir=='x'  %  parallel  case  (x  directed) 
if  diptype=- s' 

EFth=2*cos(phi).*w.*sin(beta*h*v); 

EFph=-2  *  sin(phi) .  *  sin(beta*h*  v) ; 
end 

if  diptype=- h' 
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EFth=2*cos(phi).*w.*sin(beta*h*v).*cos(pi*u/2)./(l-u.A2+le-5); 

EFph=-2*sin(phi).*sin(beta*h*v).*cos(pi*u/2); 

end 

end 

AF=sin(Nx*psix/2)./sin(psix/2).*sin(Nz*psiz/2)./sin(psiz/2); 

caftheta=abs(AF.  *EFth)/2/Nx/Nz; 

cafphi=abs(AF.*EFph)/2/Nx/Nz; 

cafthdb=20  *  log  1 0(caftheta+ 1  e-5 ) ; 

cafphdb=20*log  1 0(cafphi+ 1  e-5); 

caf=sqrt(cafphi.A2+caftheta.A2); 

cafdb=  1 0*  log  1 0(caf.  A2) ; 

figure(2),clf 

subplot(121) 

polardb(phi/rad,cafthdb,'-') 

title(['Polar  dB  plot  EYtheta,  \theta=',num2str(theta/rad),’Ao’]) 
subplot(122) 

polardb(phi/rad,cafphdb,'-')%  plot  for  E  phi  figure  2 
title([’Polar  dB  plot  E_\phi,  \theta=',num2str(theta/rad),'Ao’]) 
hold  off 

%  clear  theta  phi  AF  caf  cafdb  EFth  EFph  psi 
%  3-D  plot  % 

iplt3d=l;  %  skip  3d  plots  if =1 

if  iplt3d=0 

it=0; 

for  th=0:3:180 
jt=0;  it=it+  l ; 
for  ph=0:3:360 
jt=jt+l; 
thr=th*rad; 
phr=ph*rad; 
phd(it,jt)=ph; 
thd(it,jt)=th; 
ct=cos(thr); 
st=sin(thr); 
u=sin(thr)*cos(phr); 
v=sin(thr)  *  sin(phr) ; 
w=cos(thr); 

xsiz=beta*dz*cos(thr)-alpha; 

xsix=beta*dx*u; 

E=0; 

if  v>=0;  %  upper  half  space  (lower  half  space  zero) 

E=sin(Nz*xsiz/2)/sin(xsiz/2)*sin(Nx*xsix/2)/sin(xsix/2); 
if  E==NaN,  E=Nx*Nz;  end 
end 

if  dipdir=='z' 

EFth=2*sqrt(l-wA2)*sin(beta*h*v);  %  colinear  dipoles  (z  directed) 
EFph=0; 
end 

if  dipdir=='x' 

EFth=2*cos(phr).*w*sin(beta*h*v);  %  parallel  case  (x  directed) 
EFph=-2  *  sin(phr)  *  sin(beta*h*  v) ; 
end 

ffth=abs(EFth*E); 

ffph=abs(EFph*E); 

fth(it,jt)=ffth; 
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fph(it,jt)=ffph; 

Xth(it,j  t)=fth(it  ,j  t)  *u; 

Y  th(it,j  t)=fth(it,j  t)  *  v; 

Zth(it,jt)=fth(it,jt)*w; 

Xph(it  j  t)=fph(it  ,j  t)  *u ; 

Yph(it,jt)=fph(it,jt)*v; 

Zph(it,j  t)=fph(it  ,j  t)  *  w ; 
end 
end 

%plot  theta  component 

MAX=max([max(max(Xth)),max(max(Yth)),max(max(Zth)),max(max(Xph)),max(max(Yph)),ma 

x(max(Zph))]); 

figure(3),clf 

mesh(Xth/M  AX,  Y  th/M  AX,Zth/M  AX) 

axis([-l,l,-l,l,-l,l]) 

axis  square 

hold  on 

%  spread  N  elements  along  the  z  axis  between  [-1,1] 
z0=  [-(Nz- 1  )/2 :  (Nz- 1  )/2]  * .  75/Nz; 
xO=[-(Nx- 1  )/2 :  (Nx- 1  )/2]  *  .25/Nx; 

[X,Z]=meshgrid(xO,zO);  Y=ones(size(X))*.05; 

plot3(X,Y,Z,'ko’) 

xlabel('X=|E|sin(t)cos(p)') 

ylabel('Y=|E|sin(t)sin(p)') 

zlabel('Z=|E|cos(t)’); 

if  dipdir=='x’,  title('E_\theta  for  x  directed  dipoles  (parallel  case)’),  end 
if  dipdir=- z’,  title(’E_\theta  for  z  directed  dipoles  (collinear  case)’),  end 

%plot  phi  component 
figure(4),clf 

mesh(Xph/M  AX,  Y  ph/M  AX,Zph/M  AX) 
axis([- 1 , 1 ,- 1 , 1 ,- 1 , 1  ]) 
axis  square 
hold  on 

%  spread  N  elements  along  the  z  axis  between  [-1,1] 
z0=  [-(Nz- 1  )/2 : (Nz- 1  )/2]  *  .75/Nz; 
xO=[-(Nx- 1  )/2 :  (Nx- 1  )/2]  *  .25/Nx; 

[X,Z]=meshgrid(xO,zO);  Y=ones(size(X))*.05; 

plot3(X,Y,Z,'ko’) 

xlabel('X=|E|sin(t)cos(p)’) 

ylabel('Y=|E|sin(t)sin(p)’) 

zlabel(’Z=|E|cos(t)’); 

if  dipdir=='x’,  title('E_\phi  for  x  directed  dipoles  (parallel  case)’),  end 
if  dipdir=- z’,  title(’E_\phi  for  z  directed  dipoles  (collinear  case)’),  end 
end 
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APPENDIX  C.  MATLAB  CODES  FOR  BILOTTFS  METHOD 


This  Matlab  code  was  used  to  compute  the  voltage  ratio  versus  frequency  using  the  formulas 
derived  in  Chapter  V. 

%  FM  quadrature  demod 
%  from  note  by  Coffing  and  Main 
clear 

rad=pi/180; 

L=.45e-9;  Cp=5e-12;  Cs=5e-12;  Rp=560000; 

cc=Cs/(Cs+Cp); 

wn=  1  /sqrt(L  *  (Cs+Cp)) ; 

Q=Rp/wn/L; 

fn=wn/2/pi; 

disp(['Cs  =  ',num2str(Cs)]) 
disp(['Cp  =  ’,num2str(Cp)]) 
disp(['L  ,num2str(L)]) 
disp(['fn  =  ',num2str(fn/le6),'  MHz']) 
disp(['Q  =  ,num2str(Q)]) 

%  frequencies  in  MHz 

fl=2365; 

f2=2380; 

df=l; 

nf=floor((f2-fl  )/df)+ 1 ; 
for  i=l:nf 

fmhz(i)=f  1  +(i- 1 )  *  df; 

f=fmhz(i)*le6; 

w=2*pi*f;  %  omega 

s=j*w; 

ww=w/wn; 

wp=w+wn; 

wm=w-wn; 

H(i)=j*Q*cc*ww/(l+j*Q/ww*(wp*wm/wnA2)); 

a=2*Q*wm/wn; 

Hl(i)=Cs*L*sA2/((Cp+Cs)*L*sA2+L*s/Rp+l);  %  equation  (2) 
end 

Hm=abs(H); 

Ha=angle(H)/rad; 

Hlm=abs(Hl); 

H 1  a=angle(H  1  )/rad; 

figure(l) 
subplot(121) 
plot(fmhz,Hm) 
xlabel('Frequency,  MHz') 
ylabel('|v_2/v_l|,  dB') 
grid 

subplot(122) 
plot(fmhz,Ha) 
xlabel('Frequency,  MHz') 
ylabel('Angle  v_2/v_l,  Degrees') 
grid 
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figure(2) 
subplot(121) 
plot(fmhz,Hlm) 
xlabel(’Frequency,  MHz’) 
ylabel('|v_2/v_l|,  dB') 
grid 

subplot(122) 
plot(fmhz,Hla) 
xlabel(’Frequency,  MHz’) 
ylabel(’ Angle  v_2/v_l,  Degrees’) 
grid 
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